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Abstract 

The Nambu-Jona-Lasinio and Sakai-Sugimoto models are juxtaposed, focusing on 
the models' dynamically generated masses and the phase diagrams. The models are 
studied in the parameter space of temperature, constant electromagnetic fields and the 
quark chemical potential. When the chemical potential is zero or small, the models 
show good agreement. However, the responses of the models under the large chemical 
potential are very different. In the presence of the background magnetic field, while the 
large chemical potential triggers the de Haas- van Alphen oscillation and the boundstate 
destabilization effect of the magnetic field in the NJL model, the Sakai-Sugimoto model 
completely fails to capture these phenomena, indicating the lack of the fermi sphere in 
the latter model. Also the orders of the chiral phase transitions are generally different. 
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1 Introduction 



It has been a while since QCD, the color dynamics, was generally accepted as the theory of 
strong interaction. Our firm confidence in it mainly comes from the remarkable successes 
at high energy. In the regime of low energy, however, the degrees of freedom are apparently 
different from the QCD fundamental particles and we still do not quite understand the 
mechanism that governs this change in the degrees of freedom. The essence of the problem is 
the strongly coupled nature of QCD at low energy and the lack of powerful and efficient tools 
to tackle such theories. Currently, the best tool available is the lattice simulation of QCD 
and it has been providing affirmative results for QCD even in the low energy regime. Though 
the lattice simulation is powerful, it has its own weaknesses. For example, it suffers from the 
notorious "sign problem" when chemical potentials are included in the theory. Therefore, 
different approaches to the low energy QCD are both essential and important. One such 
approach is the construction of low energy effective theories, guided by the symmetries of 
QCD. The famous models along this line are the Nambu-Jona-Lasinio (NJL) model [1] 
and the chiral perturbation theory, and they are both known to be phenomenologically 
viable. The prominent feature of the NJL model, in particular, is the dynamical breaking 
of the chiral symmetry. Because of this feature and the sheer simplicity, this model has 
been popular, despite the model's deficiencies, such as the lack of the color gauge fields 
and confinement, also the model is nonrenormalizable and the cutoff must be explicitly 
introduced part of the definition of the theory. 

On the other hand, a very different method to handle strongly coupled field theories 
emerged as the AdS/CFT correspondence or holography pKHJIlj. It is the duality between 
string theory on the AdS space and a conformal field theory. In the low energy limit, the 
duality boils down to the equivalence between a gravitational system and a strongly cou- 
pled conformal field theory, and therefore, to obtain the information about the latter, one 
can instead study the gravitational theory, which can be easier to analyze. The dual field 
theory in the original correspondence is not like QCD, and the search for the holographic 
dual of QCD has been naturally conducted intensively. One of the most attractive models 
out of the effort is the Sakai-Sugimoto (SS) model [5]. What makes this model attractive 
is the intuitive geometric implementation of the dynamical chiral symmetry breaking. Fur- 
thermore, it is remarkable that one can obtain the spectrum of the theory, such as mesons, 
in relative ease. Therefore, even though the model does not exactly correspond to QCD 
in many parameters, such as the number of the color degrees of freedom, the spectrum 
of the SS model has been naively compared with the QCD experimental and lattice data. 
Surprisingly, the data agree fairly well, and it is currently popular to suppose a sort of 
"universality class" among QCD-like strongly coupled theories with various differences in 
parameters. 

Even with those successes, it is difficult to directly compare the SS model and QCD 
in more details, such as the response of the theories under external electromagnetic field 
and chemical potentials. This, of course, is because the strongly coupled QCD cannot 
be solved analytically and the corresponding experimental setup is hard to achieve. This 
makes it hard to tell to what extent the holographic model captures the features of QCD. 
In other words, it is unclear how universal the universality is. Therefore, in this paper, 
we turn to the NJL model, instead of QCD, to contrast with the SS model. The main 
motivation is the dynamical chiral symmetry breaking of the models and this allows us 
to investigate the chiral phase diagrams of the theories in various external parameters and 
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compare. Additionally, one can compute the dynamically generated masses in both theories 
and collate the response of the masses under different external conditions. We focus on 
those two properties of the models in the comparison. In order to dig out the similarities 
and dissimilarities, we scrutinize the responses of the models in the parameter space of 
temperature, background constant electromagnetic fields and the quark chemical potential. 
Each model has an extra parameter, the four-fermion coupling for NJL and the asymptotic 
inter-brane distance for SS, and we examine the models with respect to those parameters 
as well. 

We consider both theories for Nf = 1, this is because of its simplicity, and more im- 
portantly because the constant background fields introduced in the SS model in this case 
exactly correspond to the electromagnetic fields (we will mention this point more carefully 
in Section [3]). However, we remark that the generalization to Nf > 1 is straightforward. We 
take the large Nc limit of the NJL model and the cutoff of this nonrenormalizable theory is 
left as a parameter and we measure all the quantities in terms of this scale. Similarly, the 
quantities in the SS model are all measured in the units of the curvature scale, R. As for the 
NJL coupling, we investigate its relation to the inter-brane distance of the SS model. (The 
parameters of the NJL model are usually adjusted to be consistent with the experimental 
data. Most commonly, the NJL coupling and the cutoff are tuned to reproduce the scalar 
density and the pion decay constant from the experiment. See e.g., References [Si [7].) 

The general approaches to the models are the followings. The effective potential of the 
NJL model is derived by utilizing the worldline formalism [HI [9]. This is known as "the 
string inspired method" and the reader might feel more comfortable than other methods. 
Finite temperature is introduced by the imaginary time formalism. We pay attention to 
the subtlety with the gauge invariance in the presence of the topologically non-trivial time 
circle, and show that the careful treatment results in the automatic introduction of the 
chemical potential. For the SS model, the background fields and the chemical potential are 
introduced from the constant NSNS-i? field, similar to Reference [10]. We touch on the 
issue of the topological constraints on the constant fields and argue that the violation of 
the constraint directly leads to the instability of the system. 

Since both models are studied extensively, it is worthwhile to summarize what are known 
and what are the novel aspects of this work. First for the NJL part. Among many others, 
the work by Shovkovy in Reference [llj comes closest to our goal. (See Reference |13j for 
the similar work.) The worldline formalism is adopted in the reference for the derivation of 
the one-loop effective potential at finite temperature in the context of QED. The one-loop 
computation is exact in the large Nc limit and the resulting effective potential is identical 
to the NJL model. However, as pointed out by Gies [12], Shovkovy glosses over the subtlety 
of the gauge invariance at finite temperature. As mentioned above, we treat this issue 
carefully and reproduce Gies' result in the worldline formalism. In addition to reproducing 
the result of Reference [12], we apply the Poisson (re)summation formula to obtain the 
effective potential in more physically transparent and numerically useful form and explain 
how to treat large values of the chemical potential. Then we proceed to actually solve the 
theory numerically as the NJL model (and not as QED). 

The notable references for the worldline formalism at zero temperature are [H |l3] and 
the review paper by Schubert [9]. We mention that the more popular approach than the 
worldline formalism is the propertime formalism and it originates from Schwinger's paper 
|15j . The notable and relevant works based on the propertime formalism include Refer- 
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ences [Ml [HI [ig [H [H] . 

As for the numerical evaluation of the NJL model, the work by Inagaki et al. in Refer- 
ence |20j explores the model in the parameter space of temperature, constant background 
magnetic field and the chemical potential. Though their derivation of the effective potential 
is based on the propertime formalism, the results significantly overlap with ours. However, 
since we are not comparing the model with the experimental data but with the SS model, 
the parameter space explored in our work includes the coupling constant. Moreover, we pro- 
vide considerably deeper investigation and emphasize the importance of the fermi sphere, 
which causes the major differences between the NJL and SS models. 

Another relevant work is Reference [21] of Klevansky and Lemmer, which deals specif- 
ically with the pure electric field background at zero temperature, and the result overlaps 
with some of ours. They find the second order chiral phase transition with respect to the 
external electric field and we reproduce this result. However, in addition to the second 
order phase transition, we uncover the first order transition at a higher value of the NJL 
coupling and we believe that this first order chiral phase transition is reported for the first 
time. Furthermore, we examine the real and imaginary parts of the effective potential and 
discuss the validity of the conclusion drawn from the real part. (The imaginary part is 
not discussed in Reference [21].) We then turn on the temperature, in the background of 
the electric field, and handle this case by the weak background field expansion. Though 
we are not able to go too far in this direction, we still obtain some information about the 
dynamical mass and the phase structure. The attempt here is largely new, though the weak 
field expansion has been discussed briefly in Reference |18j . 

In above, we have listed papers of the direct relevance to our aim, but there are, of 
course, many more papers on the NJL model. For those papers, we refer to the review 
papers [Gj [7] . 

Now for the SS part. Aharony et al. have introduced finite temperature in the SS model 
in the paper [22j and we adopt their setup. The most important result of the paper for 
us is the existence of the chiral symmetry breaking in the deconfined phase. We will see 
that the NJL model has direct counterparts to the chiral behaviors of the SS model in the 
deconfined phase. The constituent mass of the holographic theories was first discussed by 
Casero et al. in Reference |23j for slightly different model from the SS model but exactly 
the same idea applies to SS and this is pointed out and examined in References \12\ 124] . 

The references listed thus far deal with the single external parameter, temperature. The 
introduction of the chemical potential in the SS model was first discussed by Horigome and 
Tanii in Reference [25] . Even though the way we introduce the chemical potential is slightly 
different from theirs, we essentially follow and reproduce their work as a special case. 

The SS model under the infiuence of the constant background electromagnetic field 
has been extensively studied by the group in Haifa [261 EH [28] and by others [291 1301 [3T] . 
Reference |26j obtains the E-T and B-T chiral phase diagrams and Reference [29] works out 
the B-T phase diagram and the response of the dynamical quark mass (m*) with respect 
to the B field. In addition to reproducing these results, we also produce ^i-T, B-fi chiral 
phase diagrams and T-m*, E-m* graphs. In Other words, we work out all the properties 
that have counterpart in the NJL model. 

The juxtaposition of the NJL and SS models is the aim of this paper and this comparative 
study is the main novelty of this work. 

Finally, the table of contents is a good guidance to the structure of this article. 
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2 NJL Model 



We start with the NJL model. Good review papers on this model are References [6l [7], 
but the following discussions will be self-contained for our purposes. As mentioned in the 
introduction, we consider the NJL model with one flavor in large Nc limit and some parts 
of the discussion are specialized for this case. We also note that throughout this work we 
use the Euclidean sign convention, including the 7 matrices. 
The NJL Lagrangian is 



,2 



^ = ^^^^ + ^§J2 + (^^75^)'] , (1) 

where we introduced the dimensionless coupling g and an arbitrary momentum scale A, 
which we will take to be the cutoff scale of this nonrenormalizable theory^ We have abbre- 
viated the summation over the color indices and explicitly, for instance, we have 

y X 2 

^ a=l ^ 

We emphasize that g is the (effective) four-fermion coupling and it does not correspond to 
the usual "strong coupling" of AdS/CFT. Comparing to the standard NJL Lagrangian, our 
Lagrangian lacks the flavor structure, because we are concentrating on the one-flavor case. 

This Lagrangian has the U{1)l x U{1)r symmetry, that is, it is invariant under the 
transformations 

qL,R ^ e''>^-^/^qL,R , qL,R ^ e-''^^-'^/^qL,R , (3) 
where are arbitrary real numbers and we have defined 

gi,/? := ^(1 =F75)9 , := 9^(1 ±75) , (4) 

with the convention 7I = 1. With those definitions, we have 

qq = qRqi + qiqR , qii^q = -iimqL - qiqR) , (5) 

and above relations make it clear that the invariance is achieved only when we have both 
quartic interactions of Equation ([T]); each quartic interaction is not chirally invariant. 

For the special case where we have Oy '■= Ol = Or, the symmetry is the diagonal vector 
part U{l)v and the corresponding transformation is 

q ^ e ^' q , q ^ qe ^' , (6) 

which leads to the quark number conservation. On the other hand, when 9a '■= —Ol = Or, 
we have the U{1)a transformation 



^ Commonly, the coupling of the four-quark interaction is defined to be G = g^K~^ /2. We find it more 
convenient to introduce the dimensionless coupling g with the cutoff scale. We will adopt the worldline 
formalism, so A will be the "propertime cutofi^" scale, as opposed to the three- or four- momentum cutoff, 
and this choice is a part of the definition of the theory. 
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This invar iance, if exists, gives the axial current conservation. 

We now want to introduce the background abehan gauge field and its interaction with 
the quarks, with respect to the U{l)v symmetry. This means that we include the kinetic 
term of the gauge field in Equation ([1]), though it is not dynamical, and replace the ordinary 
derivative with the covariant derivative 

D^ = df, + iA^ , (8) 

where we absorbed the coupling constant in the gauge field. We then have the Lagrangian 

£ = q^I^q + \{F,,f + [{qqf + {qij.qf] . (9) 

Note that the sign of the -F^ term is for the Euclidean signature which we are adopting. The 
explicit expression of F is given in Appendix [B] and we investigate the cases with constant 
electric and magnetic fields. 

It will be convenient to introduce the auxiliary fields a and vr so that the Lagrangian 
appears as 



£ =q{ilj) + ga + g^i-i^)q + \{F^u? - \N,A\a'' + vr^) . (10) 



Upon path- integration with respect to the fields a and vr, we recover the original Lagrangian 
and we have the relations 

" iv5v2^^^^ "^^^ " iV^^^'^^^^ ■ ^^^^ 

2.1 Axial Anomaly 

The key feature of the NJL model is the dynamical breakdown of the chiral symmetry, and 
for our simplified one-flavor model, the corresponding symmetry should be U However, 
since we have included the gauge fields, this symmetry does not exist in generalU In fact, 
the divergence of the axial current is proportional to FF, where F is the dual field strength 
tensor as defined explicitly in Appendix [Bj Therefore, our model is meaningful only when 

FF oc E ■ B = , (12) 

where E and B are the background electric and magnetic fields, respectively. 

When we carry out the actual numerical evaluation, we concentrate on this case. How- 
ever, up until then, we discuss the model in generality, because it is easy to generalize to 
the case with Nj > 1 and for the sake of the comparison with the SS model. 

To see what would be the order parameter of the U{1)a breaking, we note that the 
bilinears qq and qi^^q transform under ([7]) as 

qq — > cos OaQQ — sin OAqi^^q and qij^q — > sin 9Aqq + cos OAqi^^q ■ (13) 

We thus see that both a and vr are good order parameters, but what we will do is to use the 
symmetry to rotate away vr and consider (the expectation value of) a as the order parameter 
of the dynamical U{1)a breaking. 



The large does not cure this problem, because the gauge fields in focus are not the color gauge fields. 
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2.2 Effective Action 



Given a theory, all the physical information can be conveniently extracted from its effective 
action and this is especially powerful in determining the groundstate of the theory. We 
derive the effective action for our NJL model up to the one- loop order, which is exact in 
the large limit. There are several ways to construct the effective action, but we closely 
follow the classic paper by Coleman and Weinberg [32] (see also Reference [16]). Provided 
that we avoid the anomalous case, the effective action, of course, is invariant under the full 
U{\)l X U(1)r symmetry. That is, it depends only on ((t^ + tt^), therefore, as mentioned at 
the end of the previous subsection, we set vr = and derive the effective action with respect 
to the classical value of a. 

We first quickly review how the effective action can determine the groundstate of the 
theory. Define the connected (Euclidean) Green's functional W[J] by the relation 



~W[J] 



j VqVqVaexp[- j dx^{C + Ja)] , (14) 



where C is given in Equation (llOp . without the vr field, and the function J{x) is the source 
for the field a{x). Since we are interested in a constant background electromagnetic field, 
we have omitted the path integral over A^. We then have 

oJyx) dJ[x) 

where is the groundstate in the presence of the source J. This implies that we have 
the classical field (Tci(x) 

dJ[x) 

We can think of this cJci as the conjugate of J. Hence we change the variable of the connected 
Green's function by the functional Legendre transformation 

TKi] =W[J] - I dy^a,,{y)J{y) 

= - In^ VqVqVaexp[- J dx'^{C + J{a - ad)}] . (17) 

This functional, T, satisfies 



daci{x) 6aci{x) J daci{x) 



,4 5J(y) <5Ty[J] f^i6J{y) 



(5(Jci(x) 6J{y) J 6aci{x)^''^^^'' '^^^^ 
= -J{x). (18) 

Therefore, as we turn off the source J, we obtain 

d(Tcl(x) 
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and this implies that the groundstate of the theory can be obtained by extremizing the 
functional F with respect to the classical field ad- This can, in turn, be thought to determine 
the classical field ad in the groundstate. 

Now our task is to compute the effective action (jl7p and we do this to the one-loop 
order. Note that the path integral over the quark fields are quadratic and can be carried 
out exactly, yielding the factor 

In det(i^ + ga) = ^ In det[{I^f + ^V^] , (20) 

where we used the usual "75-trick". Now, since we have j^^u = ^^lu + [linlu]/'^, we get 

(^)2 =D^ + ^[j^,^^]D^D, = D' + ^j^^,[D,„D,] 

=D^ + '-7f.7uF^u , (21) 

where we used the fact that [D^, D^] = iF^^. Using this relation to Equation (pOj) . we end 
up with the expression for the effective action 

T[ad{x)] =- In jvaexp - j dx^\^^{F^u? - ^iVcAV' 



+ ^ tr ln[D2 + ^7^7^F^^ + yV^] + J{a - aci)| 



(22) 



where tr is over the Clifford algebra representation space, the color space and the spacetime 
momentum. Notice that the extra minus sign from the fermion loop makes the sign in front 
of the "trace- log" plus. 

Now as Nc becomes large, the path integral above is more dominated by the saddle 
point value at 

(J = (Jc\ ■ (23) 

This statement becomes exact at N^. = 00, on which case we are focusing. In general, we 
need to perform the derivative expansion at this stage with respect to (Td{x)^ as is done in 
the paper by Coleman and Weinberg [32] . However, since we are interested in the dynamical 
quark mass, we only need the leading order in the expansion, i.e., the constant Cci- In other 
words, we only look for the translationally invariant groundstate. The effective action F 
now is a function with respect to the c-number cjci and to this end, the spacetime integral 
in the effective action becomes trivial. Thus, it is common to define the effective potential 
via the relation 



FM =-y d^xFeff(fTcl) , (24) 

yielding 

VesM = -\{F^u? + ^iVcAV^i - 1 tr ln[I)2 + '-^^^,F^, + g^al\ . (25) 



Then the functional relation (jl9p reduces to a simple differential equation 



9 




+ 



■(g) = 



Figure 1: The tadpole cancellation condition. The dotted and solid lines represent the s and the quark 
fields, respectively. The second term on the left-hand side comes from the "vertex" NcA'^acis{x). 



Another way to obtain the effective potential is to utilize the background field method, 
by shifting the field as ct(x) s{x)+ac\ and integrating out the quark field. Then it is clear 
that Equation (j26p is precisely the tadpole cancellation condition as shown in Figure [TJ In 
what follows, we refer this condition as the "tadpole condition". 

2.3 Worldline Formalism 

We would like to evaluate the last term in the effective potential (j25p and for the convenience, 
we redefine this term as 

V{m) := ^ tr ln[D^ + '--ff^j^F^,^ + m^] , (27) 

where we recall that the trace is over the color, Dirac and loop momentum spaces, and we 
have defined the mass parameter 

m := gad . (28) 

To evaluate the effective potential, we utilize the worldline formalism that goes along the 
ideas of string theory |8|. 

To begin with, we express the function V in Schwinger's propertime [15]. By using the 
relatioril 

/■oo g-Xs 

InX = - / ds , (29) 

Jo s 

we can write 

Z Jq S I 

1 /"°° ds 2 i 
= -^nJ -e-'"" nr(x|exp[-s(D2 + -^^^,F^,)]\x) , (30) 

where the operator tr now is only over the representation space of Clifford algebra and the 
factor A^c came out of the trace over the color space. 

Now in the worldline formalism, we want to regard as fields in the one-dimensional 
space of circular worldline, parametrized by, say, r, and express the effective potential in 
terms of the path integral over the fields x^(t). The proper length of the circle's circumfer- 
ence is the parameter s, which is called "modulus" in the string literature, and the factor 
1/s in V removes the redundancy from the translation and reversal of the worldline coor- 
dinate. In order to deal with the last trace over the Clifford representation space, we need 

^ This is the analytically continued version of the propertime and is the "imaginary propertime" . Later, 
we will come back to this point concerning the real and imaginary propertime. We note that this is essentially 
different from the analytic continuation of the spacetime coordinate. 
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one more ingredient. What we need to do is very similar to the treatment of worldsheet 
fermions in superstring theory. We introduce constant Grassmann fields ip^ that satisfy 

{'iPii,i^u} = Sf,u , (31) 

and introduce a basis of the Clifford representation space \a) so that 

i^^la) = -L(^^)^^|/3) . (32) 

We then can write 

V{m) = -lNc I ^e-"'^'{a,x\exp[-s{D'^ + it^f,'ilj^F^u)]\a,x) 

^e-'^'^M I^Vx{t) I^ViP{T)ew[- dT{C,+C^)] , (33) 

with 





/■°° ds 








Jo s 












Jo s 



.e' 



Cx ■■= -^x +iA- X and := -ip ■ ip - i—ip ■ F ■ ip . (34) 

Here the dots on the fields denote the derivative with respect to r and e' is the einbein 
of the worldline circle, which can be an arbitrary fixed number because we have decided 
to encode the moduli in the parameter s (see, for example, Chapter 5 of Polchinski |33j). 
Hence from now on, we set e' = 2, which is the common choice in the literature. As usual, 
the factor J\f arises from the integration over the conjugate momenta to obtain the path 
integral in the Lagrangian form, and it satisfies the relation 

N f Pxexp[- / dT-iP'] = (47rs)-'^/2 ^ (35) 
Jp Jo 4 

where d is the spacetime dimension, i.e., 4 in our case. Finally, the subscripts P and A on 
the path integral symbols imply the boundary conditions of the fields 

x{0) = x{s) and ^(0) = -^'(s) , (36) 

for all the components. 

Before we proceed further, we discuss the treatment of finite temperature in the context 
of the worldline formalism. 



2.4 Finite Temperature and Chemical Potential 

As we have adopted the Euclidean signature, we use the imaginary time formalism to handle 
finite temperature. This means that the time direction in the spacetime is compactified 
with the circumference /3, and a worldline can wind around this direction. Thus, the path 
integral should be carried out in each winding sector, and then they must be summed with 
appropriate statistical factor [34] (see also Refs. ^SSjiSGj). This can be written as 

oo 

^xx = ^ i~^TP^xo+nl3,Xi) ' (37) 
n=—OD 
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p 




Figure 2: Conceptual picture of the path integral at finite temperature. The left-most loop paths corre- 
sponds to n = 0, I.e., the T = contribution. 

where Pxy represents the path integral over the paths from a spacetime point x to y and 
poo jj^eans the path integral in M^, i.e. without the topological condition. We emphasize 
that the alternating sign (— l)" comes from the spacetime fermions q, and they should not 
be confused with the worldline fields. See Figure [2] for a conceptual picture of the path 
integral. 

Let us look at the path integral of x(r) in Equation (f33l) . We have the Lagrangian Cx 
as shown in Equation ([M|) and this does not appear to be quadratic in the field x^{t). At 
this point, it is very common that a specific gauge, called the Schwinger-Fock gauge: 



is adopted to bring the Lagrangian into the quadratic form. However, one immediately 
realizes the problem with this gauge at finite temperature by re-examining Schwinger's 
original paper [15]. There, the holonomy factor, 



appears in the two-point Green's function, while in the gauge psp this factor is absent. 
This factor has no relevance if one is interested in the effective potential where x = y, and 
if the spacetime is topologically trivial. In the imaginary time formalism of the theory at 
finite temperature, the latter condition is false. In the geometrical language, the transition 
function defined as a map from the time circle to the U{1) gauge group is non-trivial 
and classified by 7ri[?7(l)] = Z. This in turn implies that the gauge transformation is 
severely restricted to a certain class and one must be very careful with the gauge fixing. In 
other words, the consequence of the time circle is (geometrically) exactly the same as the 
charge quantization of the Dirac monopole and the Aharonov-Bohm effect (see for instance 
Reference [37]). As in those cases, the parallel transport of the wavefunction around the 
time circle, i.e., the Polyakov loop, obtains a non-trivial physical significance. As we will 
see below, this automatically introduces the chemical potential into the system. 

Gies |12] appears to be the first one to pointed out this subtlety and he generalized 
Schwinger's propertime formalism to finite temperature with the careful treatment of the 
holonomy factor. We are going to show how this works in the worldline formalism. 

To take advantage of the Schwinger-Fock gauge, we split the Lagrangian into the 
Schwinger-Fock gauge part and the rest 




(38) 




(39) 



Cx = \x^ + ix ■ {A + -F ■ X -- 
4 ^ 2 2 



F -x) 



( 



1 . 



-X ■ F ■ x\ +ix ■ {A + -F ■ x) . (40) 
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Now the important point is that 



d{A + ^F-x) = 0, (41) 



oo 



and this imphes that within a homotopy class, the integral 

/ drx ■ {A + -F ■ x) = dx-{A + -F ■ x) (42) 

Jo 2 y^(o) 2 

is independent of the path. Therefore, this factors out of the path-integral. Furthermore, 
without loss of generality, we may evaluate the x integral of Equation ()42p along the straight 
line in xo-direction. In this case, the integral over F ■ x vanishes due to the anti-symmetry 
of F and we are left with the integral of ^o- In what follows, we are going to consider the 
case with constant and we define a real number 

H:=iAo. (43) 

Then we can write 

^ Jo ~ 

xM Vx{t) / PV'(T)exp{- / dT{C^sF + C^)] , (44) 

Jx(0) J A Jo 

1 i 1 • 

^xSF ■= -^x^ - -X ■ F ■ X and := -ip ■ tp - iil^ ■ F ■ if^ . (45) 

The subscript "SF" of the bosonic Lagrangian signifies the form in the Schwinger-Fock 
gauge. The integer n is the winding number in the xq direction and it is also implicitly in 
x{s) = x(0) + nPxQ, where xq is the unit vector in the time direction. 

2.5 The Effective Potential 

In this subsection, we derive the expression for the effective potential that is ready for the 
numerical evaluation. 

In Equation (j44p . we are interested in the constant F, i.e., it is independent of the 
paths x{t). Thus we can separately calculate the path integrals for the worldline bosonic 
and fermionic parts, 

f'X^^s'j PS 

Ib :=M / Vx{t) exp ( - / drC^sp) 

Jx{0) Jo 

If:= [ V^{T)exp{- [ drC^) , (46) 



with 



where the factor naturally comes with the bosonic part because of the relation p5|) . The 
evaluation of them is rather lengthy and is given in Appendix [Al The results are 



^1/2 

Ib =(4^s)~2exp [- (n/3/2)2{Fcot(Fs)}oo] ( det^ 



sin(Fs) 
Fs 

If =4(detL[cos(Fs)])i/2 , (47) 
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where the functions of the matrix F^^ are defined through the formal power series and det^ 
imphes the determinant with respect to the matrices with the components in the Lorentz 
indices. We notice that the Lorentz invariance is exphcitly broken in Ib, due to the finite 
temperature. 

At this point, V has the form 



V{m) 



87r2 



dss 



XdetL[cos(Fs)])^/2( detL 



sin(Fs) 
Fs 



-1/2 



^ (-l)"e-"^'^exp [ - (n^/2)2{Fcot(Fs)} 



00 



(48) 



In this form, the physical interpretation of the parameter is still not clear. Moreover, it 
turns out that this form is not suitable for the numerical evaluation in a certain range of 
parameters. Thus, we use the Poisson (re)summation formula 



exp(— vran^ + 27ri6n) = a ""^^^ exp [ (n — 6)^] , 



(49) 



to convert the above expression to 

T 



V{m)/N, 



47r3/2 



dss-^/^{{Fs) cot(Fs)}oo^/^{detL [(Fs) cot(Fs)] } 



1/2 



^ exp [ - {{Fs) cot(Fs)}oo^(27rr/ - ijifs 



(50) 



^1/2 



where T := 1/(3 and I runs over the half integers '^1/2 ■= {n + 1/2 : n E Z}. In this 
form, we see that the parameter fi is the chemical potential, as advertised before. In the 
standard method, the chemical potential is included from the beginning as a part of the 
grand partition function (in the end, this amounts to the shift Po ^ Po — ifJ')- Contrast to 
this, the automatic inclusion of the chemical potential we see here is rather remarkable. 

We now have to evaluate the terms involving the matrix F. The determinant can be 
calculated relatively easily by using Schwinger's method [K], but {(Fs) cot(Fs)}oo is not 
straightforward and for this, we adopt the method of Reference [38] (see also |39|). They 
are computed in Appendix [B] and the results are 



{detL[(Fs)cot(Fs)]} 



1/2 



gcoth (^s^.F+ v^2 + g2^ cot (^s^-.F+/.F2 + g2^ 



(51) 



and 



{(Fs)cot(Fs)} 



2V.F2 + g2 

I -i(52 + F2) + V.F2+g2|^^+^^2+gicoth (^s^/7W^^+^^ 



+ |i(52 + F2) + V^2+g2|^_^+/^^7^COt(s 



(52) 
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where the Lorentz invariants are defined as 

^:=-^F^,F^, = ^(^2-^2) g _'_jp^^p^^ ^ g . g ^ (53) 

Those expressions are very clumsy, so we consider a few special cases. As discussed in 
Section [XTl our Nf = 1 model does not make sense unless G = due to the U{1)a anomaly. 
The case with ^ 7^ is of interest when Nf > 1, but in what follows, we concentrate on the 
non-anomalous case. 

At zero temperature, the Lorentz invariance should be recovered and we expect the 
expression for the effective potential becomes considerably simple. In fact, we can go back 
to Equation (|48p and simply set n = 00 We then have 

AT rOO 

V{m) = J dss'^e^"'^'{sV2Tcoth{sV2^)} . (54) 

From this, we see that for the cases with \B\ > \E\ and \B\ < \E\, the system behaves as if 
there is only i?-field and -E-field, respectively. Also as is well-known, the pure E'-field case 
can be obtained from the pure -B-field case just by sending B — > iE, and vice versa. 

When we turn on the temperature, all those features are lost. Examining Equation (I52p . 
one sees that the response of the system for the cases \B\ ^ \E\ is not as simple as the zero 
temperature case. Moreover, the pure B and E background systems are not related through 
B — > iE. These are all due to the lack of the Lorentz invariance at finite temperature. 
Because of this complexity, we find it best to handle the finite temperature system with 
only B and E fields separately. 

2.5.1 Pure B Background 

We have for this case 

{detL[{Fs)cot{Fs)]y^'^ ^ s\B\coth{s\B\) and {(Fs) cot(Fs)}oo ^ 1 , (55) 
which lead to the expression of the effective potential 

Vesim) =-B^ + -N.k^g-^m' + ^ dss->/\s\B\ coth(s|i?|)} 

X ^ eyi-p[- s{{2'kTI -ifif + rr?]] . (56) 

First, notice that the first term on the right-hand side is A^c-suppressed, therefore, we 
drop this term from now on. Next, observe that the s integral is ill-defined at s = and 
we regulate this by introducing the propertime cutoff s = 1/ for the lower bound of the 
integral0 (Up to this point, the scale A has been an arbitrary momentum scale but now it 

* At zero temperature, the dependence on the chemical potential drops out, as the time direction is 
noncompact, and introduction of the chemical potential at zero temperature requires the standard method 
mentioned previously. We will generally consider the zero temperature case without chemical potential. 

^ Reference f2T| assumes that the cutoff A can depend on the parameter m. This assumption requires 
another equation, in addition to the tadpole condition (|26|l . For this, they derive a consistency condition, 
which corresponds to Equation ((18} in our path-integral treatment of Section [2.21 However, the large Nc 
makes the classical path totally dominant in the path integration as discussed around Equation (|23|l . This 
results in the fact that the effective action is independent of the source J and the relation (llSp gives nothing 
other than the tadpole condition (|26p . still lacking an equation that determines the function A(m). It seems 
that the consistency condition of Reference [JT] comes out of the lack of the proper change of variable via 
the Legendre transformation, in order to obtain the energy of the system. 
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is taken to be the propertime cutoff scale.) We would like to measure all the quantities in 
the units of the cutoff scale A. This means that we rescale all the dimensionful parameters 
with A, for instance, s — > sA^ and B — > B / . Then defining the dimensionless quantity 
V^flf) we have 

Ves{m) — 



=^ + ^^ j dss-3/2coth(s|S|) ^ exp[-4(27rT/-i/x)2 + m2}] , (57) 

where all the parameters are now dimensionless. 

Though we have regulated the lower end of the s integral, it is finite at the upper end, 
only when 

Re [{2tiTI - ifj,f + = {2TiTlf - ^? + m'^ > Q . (58) 

This condition crucially depends on the value of the chemical potential. It is well-known 
that the chemical potential of free (complex) scalar field theory may not exceed the mass of 
the particle, because the chemical potential effectively acts as a negative mass squared. One 
might naively think that the similar situation is occurring in our theory when the chemical 
potential is too large. However, this is false. We must recall that we have been adopting 
the analytically continued "imaginary propertime" in Equation ()29p . We have just found 
that this analytic continuation is appropriate only when the condition (j58p is met. If we 
consider Equation ()29p as the continuation to the positive imaginary axis, there is nothing 
that prevents us to continue to the negative imaginary axis to make the propertime integral 
well-defined. In other words, when the condition ()58p is violated, we must analytically 
continue to negative s. 

This is the idea of treating the case with a large chemical potential, but there are 
complications due to the existence of the cutoff, the cut in the complex s-plane from the 
factor of s~^/^ and the poles from the hyperbolic cotangent along the imaginary s-axis. We 
work out the details in Appendix [Cl The result is 

%s = 

™2 T-T r /"OO 

" j ds(sS) coth(sS)s-5/2^ exp [- s{(27rr/ -i;u)2 + m2}] 



-I Re 



l>l 

I 



/oo 
ds{sB) coth(sS)s"^/2 ^ exp [s{{2t:TI - i^f + rr?]] 

1=1/2 

+ i r d^{e'^B)coi\i{e'^B)e-'i'^ V exp [ - e^'^{(27rr/ - i/i)^ + m^}] 
^0 



1=1/2 



I oo / D \ 3/2 



1=1/2 k=l 



(59) 



where B := \B\ and we have the condition B < n, indicating that B may not exceed the 
cutoff scale squared A^ too much. The parameter / is defined to be 

l■.= [{2^TTy^e{^l-m)^/J]?^^]^, (60) 
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where without loss of generahty /x and m are assumed to be greater than or equal to zero, 
6{x) is the Heaviside theta and the symbol [x\g is the half-integer version of the Gauss 
symbol, ie., it is the largest half integer less than or equal to x. Obtaining Equation ([59]) 
involves contour integrations, and referring to Figure [l2l the first line is the contribution 
from the positive horizontal axis, the second line is from the negative horizontal axis, third 
one is from the small circular contour (coming from the existence of the cutoff) and the last 
term is the contribution from the poles mentioned above. One can see that the last term is 
oscillatory and we will see that this part gives rise to the well-known de Haas-van Alphen 
effect [17j and it complicates the phase structure. 

The effective potential obtained above contains all the necessary information and is ready 
for the numerical evaluation. Regarding V^fj as a function of m^, the tadpole condition (j26p 
becomes 

^ _ dVcK _ d{m?) dV^s - 2g^cT i '^^''^ (61) 
dad dad d{vn?) ^ d{m?) ' 

where we used the definition (j28p . We denote the solution to the tadpole condition by m* . 
We see that m* = is always a solution and the solution other than this satisfy 



d(m^ 



. (62) 



When there are more than one solution, energetically preferred one must be chosen by 
consulting with V^s- In this way, we can determine the dynamical mass and the phase 
structure of the theory. 

2.5.2 Pure E Background 

We have for this case 

{detL[(Fs)cot(Fs)]}^/^ ^ s|-E| | cot(s|^|)| 

{{Fs) cot(Fs)}oo s\E\ cot(s|^|) , (63) 

which lead to the expression of the effective potential 

Kff(m) =^ + dss-^/'^^/sEcot{sE) 

X ^ exp[- s{{sEcot{sE)y\2TTTl + m'^}] , (64) 

where as in the previous subsection, Ves and all the parameters are made dimensionless 
with respect to the cutoff A and we redefined E := \E\. Because of the cotangent in 
the exponent, we have essential singularities and this makes the examination of this case 
notoriously difficult. In the following, we make two attempts: T ^ limit and weak E-&eld 
expansion. 
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Zero Temperature 

In the zero temperature limit, instead of summing over I, we should integrate Equation (|64p 
with respect to the combination Tl. This yields 

K,ff=^ + ^y^ t^s5-'{(sii;)cot(si?)}e-'"^ (65) 

As noted in the footnote near Equation (j54p . the chemical potential drops out and we 
consider the system without it. We find that only the simple poles of the cotangent are 
present. Thus, it is commonly evaluated through the convention 

oo 

cot{sE) p. V. cot{sE) + ^ V <5(s - TTTi/E) , (66) 

n=l 

where P.V. denotes the principal value. 

The principal value gives the real part of the effective potential and this has been 
evaluated in the literature. We quote the result from Reference [21] , 

Re [Feff] = ^ + ^1 Re[Q] - ge + E^im^) - ^E^E,{m^)^ , (67) 

where 

with z := im? /{2E), 

UE ■■= t dss-^e-''^''{{sE) cot{sE) - 1 + {sEf/3} , (69) 
Jo 

and the functions En{x) are the exponential integrals defined as 

^e-^* . (70) 

Using the real part of the effective potential, Klevansky and Lemmer ^21j examines the 
dynamical mass and the phase structure of the system under the stress of the constant 
electric field. (We will reproduce their results in the next subsection.) 

The second term in the replacement ()66p yields the imaginary part of the effective 
potential 

I-[^efr]=^f;f-)V^-\ (71) 

n=l ^ ' 

The imaginary part is interpreted as the pair creation rate of the quarks in Reference |15j . 
The pair creation inevitably implies the instability of the groundstate and in our system, 
this occurs at any finite value of E. Therefore, the conclusions about the mass and the 
phase structure drawn from the real part may not have any physical significance. However, 
as long as the imaginary part is small compared to the real part, we can assume that the 
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instability sets in slowly, compared to the dynamical mass m*, and we can still obtain the 
physical insights. In the numerical evaluation, therefore, we will simultaneously examine 
the real and the imaginary parts of the effective potential. 



Weak E-Field Expansion 

We attempt to circumvent the essential singularities in Equation (j64p by expanding in small 
E. This, of course, is a very dangerous thing to do. However, we extract some qualitative 
conclusions by comparing it to the zero temperature case and also we check if the higher 
order terms in the expansion give less contributions to the effective potential. 
The expansion itself is straightforward and we have 

,2 



^ / dss'^^^ {l + 2s{2TTTl-inf}exp[- s{{2TTTl-ifif + m'^}] 



+0{E^) . (72) 

One can easily go up to higher orders as desired. We must evaluate the sum and the integral 
term by term in the expansion as in Appendix O However, assuming that the parameter 
|/| is large enough to satisfy the convergence condition (|58p . one can carry out the integral 
before the sum and obtain a closed form. As discussed in the appendix, one can then 
analytically continue the integral beyond the validity range of the parameter I to the whole 
half integers and this considerably eases the practical evaluation^ Thus in practice, one 
should exchange the order of the integral and the sum in Equation (j72p . without worrying 



about the convergence condition (j58]) rl 



2.6 Numerical Evaluation 

We are interested in the dynamical mass m* and the effective potential V^fi in the five- 
dimensional parameter space (T, fi,B,E,g). It is, of course, not practical to cover entire 
space at once, so we examine slices of the parameter space. As a warm up and to illustrate 
the general behavior of the NJL model, we first investigate the slice oi B = = E. We 
then turn on those parameters one at a time. 



2.6.1 General Observation without Electromagnetic Fields 

We learn in this section that the temperature and the chemical potential both behave as 
the "destabilizer" of the quark boundstate, and that when the condition (ISSh is violated 
and the analytic continuation is required, the effective potential shows more complicated 
behavior than when the condition is satisfied. Those are the generic properties of the model 
and remain essentially the same when we turn on the external electromagnetic fields. 

® For the pure B background, we have not been able to obtain the closed form of the s integral, so needed 
to evaluate the very complicated expression in Equation (|59p . 

^ Also it is practically convenient to realize that the negative I is exactly the complex conjugate of the 
positive counterpart. 
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The tadpole condition (j26p for this case can be written as 



TT' 



3/2 



TRe 



/oo 



(73) 



We have indicated that the integral should be done before the sum. This is because the 
integral for this case is simple enough to yield a closed form, and as discussed in Appendix [Cl 
and Section 12.5.21 we can then analytically continue the integral to the values of I that 
violates the condition ([58]) . 

We first plot the right-hand side of the tadpole condition ([73]) in Figure El These plots 



fi-O fixed 



0.14 
0.12 
0.10 
0.08 
0.06 
0.04 
0.02 
0.00 




0.12 
0.10 
0.08 

CO 

0.06 
0.04 
0.02 
0.00 





= 




f 


= o3v 




f 


= 6 











0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 3: The right-hand side of Equation fTS]) . i.e., 7r^''^/g^, against m*. For T = = the value of 
RHS at m* = is l/47r^''^ « 0.141 and this corresponds to the famous critical coupling gc = 2%. Note the 
different behavior appearing for {TvTf - ^2 < 0, which is the case for T = 0.1 and ^ = 0.6. In this case, 
there can be two solutions with m* 7^ 0. 

show that there is a critical value of the coupling g, below which the solution to the equation 
does not exist (and the only solution to the tadpole condition is m* = 0), and the system 
is in the chirally symmetric phase. The critical coupling is plotted against the temperature 
in Figure [H In particular, the critical coupling for T = = is the famous gc = 27r 




0.15 
T 



Figure 4: The critical coupling gc is defined to be the value, below which the solution to the tadpole 
condition does not exist. Note that for small T and large fi, g = gc does not imply m* — 0, as shown in the 
right panel of Figure [3] The plots show that the larger T and fi require the larger coupling. 
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(^3/2/^2 ^ 0.141) of the original work We observe that both the temperature and 
chemical potential tend to destroy the quark boundstate, as these parameters necessitate 
stronger coupling. 

In the right panel of Figure [3l we notice the difference between the low and high values 
of the chemical potential. For the parameters given in the diagram, we have (vrr)^ - fi^ > 
for {T,fi) = (0.1,0.3) and negative for (0.1,0.6). Thus, when the chemical potential is 
large enough (and the temperature is low enough) to require the analytic continuation, 
the system respond in a more complicated way and in particular, there can be multiple 
solutions to the tadpole condition for a given value of the coupling. One manifestation 
of such complication is the order of the chiral phase transition. In Figure [5l the effective 
potentials are plotted against the parameter m for the low T, high fi region and for the high 
T, low fi. The reference point of the potential is taken to be the value at m = 0, which is 
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Figure 5: The plot of Al^ff Vcs{m) — Veff(O). The left panel shows the first order phase transition with 
respect to the chemical potential in the low T, high fi region and the right panel shows the second order 
transition in the high T, low fi region. 

always a minimum, and this base value defines the difference potential AT4ff. Notice that 
the low T, high /i diagram on the left has a large range of m where (vrT)^ —fi'^-\-m? < 0, and 
this leads to the extra structure, the potential barrier, that causes the first order chiral phase 
transition. Meanwhile, the high T, low /i diagram of the right panel has (vrT)^ — /i^ > for 
all the values of /i plotted. For this case, the analytic continuation does not occur and Ves 
is relatively simple, leading to the second order phase transition. The complication due to 
the analytic continuation discussed briefly here is the generic phenomenon and we will see 
the further consequences in the next subsection. 

Finally in the absence of the electromagnetic background, we show the response of m* 
with respect to the coupling g in Figure EJ As expected, the dynamical mass generally 
increases with the coupling. Diagram (b) zooms into the near critical coupling Qc, and 
shows the relation m* ~ in this region. This response of the dynamical mass with 
respect to g will be contrasted to the dynamical mass of the SS model with respect to the 
probe brane separation L. 



2.6.2 Pure B Background 

We now turn on the B field (and E = 0) and the effective potential for this case is given in 
Equation ([59]) . 
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Fi glirc 6: These are plotted for T = 0.1 and /i = 0. Diagram (b) is m*^ against the region of g near the 
critical value and shows the relation m* ~ y^. The linear behavior of (b) fails in the larger values of g. The 
graph for T = is nearly identical, except that the critical value is gc — 2tt and the graph is slightly shifted 
to the left. 



The salient feature of this effective potential is the last term in the equation, the os- 
cillation term. When the condition (j58p is violated, this term comes in effect, leading to 
the oscillatory behavior of the effective potential with respect to the various parameters. 

Figure [7| is the plot of the effective potential with respect to the parameter B and the 
other parameters fixed. We clearly see the oscillation, known as the de Haas-van Alphen 
effect, and this is caused by the crossings of the fermi sphere through the Landau levels. 
This effect leads to the multiple solutions to the tadpole condition ([26]) and complicates the 
model in the low T, high fj, region of the parameter space. 

From Equation (I59jl . it is easy to work out the tadpole condition ()26[) and it can be 
written in the form vr'^/^/g^ = RHS, similar to Equation (j73p . The right-hand side (RHS) is 
plotted for various parameters in Figure [HI The plots of the left panel show the behavior of 
RHS similar to Figure [3] without the B field. One important difference, though, is the case 
with T = where the presence of the B field makes the critical coupling 0, i.e., the magnetic 
field stabilizes the quark boundstate and any finite coupling breaks the chiral symmetry. 
For this reason, the magnetic field is sometimes called the "catalysis" of the chiral symmetry 
breaking |19j . However, even with the B field, any finite temperature brings the critical 
coupling nonzero. Figure [9] shows the plots of critical coupling against the magnetic field 
at // = 0. In the figure, we see the competing effects of B as the stabilizer and T as the 
destabilizer of the quark boundstate. 

Now, the right panel of Figure [8] is plotted for the parameters that satisfy (vrT)^ — /i^ + 
m*^ < for m* < 0.5. We can observe the effect of the oscillation, and for a certain range 
of the coupling and B {e.g., B = 0.1 and vr^/^/g^ ~ 0.07 at g'^ = 80), there are four m* ^ 
solutions to the tadpole condition. We must consult with the effective potential to decide 
which solution is the stable and energetically preferred one. Figure [10] is the plot of the 
potential. We indeed observe four peaks and troughs. The top plot is at the critical 
value of the chemical potential and above this value, the mass sharply drops to zero and 
the system goes through the chiral phase transition to the symmetric phase. Curiously, the 
bottom two curves show that there is another first order phase transition that is different 
from the chiral phase transition. During this transition, the mass sharply drops to a lower 
value but not to zero. This is the extra structure introduced by the oscillation. 
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Figure 7: The plot of the effective potential with respect to B, illustrating the de Haas-van Alphen effect. 
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Figure 8: The right-hand side (RHS) of the tadpole condition similar to Equation (|73|l is plotted as a 
function of m*. On the left panel with T = 0, the graph blows up to infinity as m* 0, indicating Qc — 0. 
Finite T brings Qc finite. The right panel is plotted for the parameters that satisfy (ttT)^ — fi^ + m*^ < 
for m* < 0.5. Notice the effect of the oscillation that leads to multiple solutions to the tadpole condition for 
a certain range of g. 



fi-Q fixed 




1 p - 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

B 



Figure 9: The plots show that the B field is the stabilizer of the boundstate, as it requires smaller value 
of the coupling at high B, and as discussed previously, the temperature is the destabilizer. 



23 



j-=80, T=0.005, B=0.1 
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Figure 10: The plot of the effective potential for nonzero B-field and in the low T, high n region of 
parameters. In addition to the potential barrier that causes the first order chiral phase transition, we have 
extra structure due to the de Haas- van Alphen oscillation that causes another first order transition, different 
from the chiral phase transition (the bottom two plots). 




Figure 11: The plot of the dynamical mass against the chemical potential /i. The left panel shows the region 
of low T, high fj., where the mass goes to zero abruptly, corresponding to the first order phase transition. 
For the nonzero B-field, we have an additional structure due to the de Haas-van Alphen oscillation and the 
system goes through a sequence of first order phase transitions. The right panel shows the region of high T, 
low n, where the mass goes to zero smoothly as the chemical potential is increased. This corresponds to the 
second order chiral phase transition. 



Figure [TT] of the /x-m* graphs clearly summarize the discussions made so far. On the 
left, the graph in the sohd hne directly corresponds to the behavior of the effective potential 
shown in Figure [10] and reveals the existence of the two first order phase transitions. In 
the same diagram, the dashed graph is plotted for B = 0, indicating that there is only one 
first order chiral phase transition (see the left diagram in Figured]). Therefore, the extra 
first order phase transition is visibly due to the de Haas-van Alphen effect at S / oi On 
the right panel, we have plotted the //-m* graph for the higher temperature so that the 
condition (|58p is always satisfied. For this case, the mass smoothly goes down to zero and 
this is the second order chiral phase transition, similar to the one shown in the right panel 

* We remark that the sharp first order phase transitions shown here is slightly different from rather 
smooth transitions in Reference [20] , The smoothness is most likely due to the higher temperature in their 
setup. 
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of Figure [5j Notice also the role of B and T as the stabilizer and destabilizer, respectively. 

We now discuss the response of the dynamical mass with respect to the B field. First, 
Figure [l2] shows the simpler case with /i = 0. For this value of the chemical potential, 
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Figure 12: The plot of the dynamical mass with respect to B a,t fj, — fixed. The temperature here is low 
enough so that the chiral symmetry is broken at i? = 0. The right panel shows the relation m* ~ for 
the region of small B. The parameter mo is defined as mo := m*{B — 0). 

the mass responds intuitively; it increases monotonically as the external magnetic field (the 
stabilizer) increases. Notice that for this plot, we have chosen the low enough temperature 
so that the chiral symmetry is broken even at i? = 0. The right panel of the figure shows 
that the mass increases quadratically with respect to B in the region of small B. Also the 
diagram shows that the tangent of the graphs increases as the temperature is raised. 
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Figure 13: The plot of the dynamical masses with respect to B. For (T, /i) = (0.3,0.3) and (0.39,0), the 
masses continuously rise from zero as B increases, corresponding to the second order phase transitions. For 
{T,ii) — (0.2,0.4), the mass goes smoothly goes to zero in the region of small B and rises abruptly in the 
region of large B. They correspond to the second and first order phase transitions (see the left panel of 
Figure [14]). We note that for a certain range of parameters, the left part of the dotted graph can exhibit a 
sequence of first order phase transitions similar to the one shown in m*-/i diagram of Figure [TT] 

The B-m* graph with nonzero chemical potential is shown in Figure [T3l Let us first 
observe the dotted curves for T = 0.2 and /x = 0.4. The response of the mass against B 
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Figure 14: The left panel shows the behavior of the effective potential for {T,fi) = (0.2,0.4), exhibiting 
the first and second order phase transitions (see Figure [TS]). The right panel is plotted for the high enough 
temperatures, concentrating on the values near B — Bo of the second order phase transition points. The 
plots reveal the relation m* ~ near B = Bo and the coefficient increases as temperature goes up. 



is highly counter-intuitive for fower B: the mass decreases as the magnetic field is dialed 
up. Moreover, this curve in the low B region can exhibit a sequence of first order phase 
transitions that are different from the chiral transition, due to the de Haas-van Alphen 
oscillation. (We have omitted this very complicated graph.) The further increase in the B 
field leads to the second order chiral phase transition around B ~ 0.27. But approximately 
at 0.78, the system comes back to the broken phase through the first order transition and 
the mass increases as B is raised. The left panel of Figure [E] shows the series of action in 
terms of the effective potential. 

The complicated behavior here requires a digression. As argued in Reference [19j . the 
stabilizing effect of the B field directly comes from the dominance of the lowest Landau level 
(LLL). In Figure [T2I we have plotted for /i = and the system is always in the LLL, leading 
to the stabilization effect. When fi is large and T is relatively small, we have a reasonably 
well-defined fermi sphere (smeared due to the temperature) and the fermi energy can be 
much larger than the Landau level spacing. If this is the situation of the system, the effect 
from the LLL is not significant and we do not expect the stabilizing effect of the magnetic 
field. Even though the finite temperature and the interaction effects prevent the precise 
estimate here, this in fact is roughly the case for the low B region of the dotted curve in 
Figure [T3l Here we observe that the B field acts as a destabilizer of the quark boundstate. 
As the magnetic field becomes large enough, the Landau level spacing becomes also large 
and the stabilizing effect from the LLL takes place. To summarize the digression, we can 
conclude that the oscillation effect and the destabilization effect of B are in general caused 
by the existence of the fermi sphere and its relation to the Landau levels. 

Let us now turn to the other plots in Figure [T3l These are the cases where the tem- 
peratures are high and the chemical potentials are low enough so that the system is in the 
symmetric phase in the region of low B. The masses rise continuously from zero at critical 
values Bqs, indicating the second order chiral phase transition, and increases monotonically 
with B. The right panel of Figure [E] shows that the mass rises from zero as y/B and the 
coefficient increases as the temperature. 

The mass depends also on the temperature and the T-m* diagram is shown in Figure [T5l 
Here we see that as long as = 0, or when it is small enough, the mass continuously goes 
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Figure 15: The plot of the dynamical mass as a function of the temperature. We see that the larger values 
of B lead to the larger mass, i.e., B is the stabilizer. For = 0, the mass smoothly becomes zero as the 
temperature is increased, indicating the second order chiral phase transition. Meanwhile, for ^ = 0.4, the 
mass drops abruptly to zero at a critical temperature, indicating the first order phase transition. 

to zero as the temperature is increased, indicating the second order chiral phase transition. 
We can also observe that the larger the B field is, the larger the mass is, i.e., B acts as the 
stabilizer of the quark boundstate. When the chemical potential is high, such as = 0.4, 
the first order phase transition occurs, as shown as the dashed curve in the diagram. 

We now proceed to discuss the phase diagrams. Though we have seen that there are 
more than one kind of phase transition, here, we focus on the chiral phase diagrams. 

The B-T phase diagram is shown in Figure [161 The top two curves for /.i = and 0.3 
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Figure 16: The B-T chiral phase diagram. The upper temperature regions are the chiral symmetric phase. 
The dotted curve shown on the left panel indicates the first order phase transition and the solid ones are 
second order. The right panel plots the small B region against {\T /Tq — 1\)^^'^ where To is the phase transition 
temperature at B = 0, showing the relation T ~ ±5^. 

shape in expected way as T being the destabilizer and B as the stabilizer. However, the small 
B region of the bottom curve for the large chemical potential jj, = 0.4 is rather surprising. 
The critical temperature decreases as B is increased, indicating that the magnetic field is 
acting as the destabilizer in that region of B. When B is further dialed up, the first order 
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Figure 18: The fi-T phase diagram. In the low temperature region, the order of the transition becomes 
first order. It is counter-intuitive that the critical value of the chemical potential in the low temperature 
region is smaller for large value of B, and this is due to the suppression of the LLL effect. 



transition is triggered and then B starts acting as the stabiHzer. We have seen the similar 
phenomenon when we discussed the B-m* diagram of Figure [T3] and noted that this is 
caused by the existence of the fermi sphere and the smaU Landau level spacing. This effect, 
in the B-T diagram, appears as the "dent" in the chiral phase transition curve. The right 
diagram of Figure [T6l zooms into the region of small B and shows the relation T ~ ±B^. 
Note that the slops of the plots decrease as the chemical potential increases. 

The "dent" observed in the B-T phase diagram can also be observed in the B-fj, phase 
diagram of Figure [TTl This phase diagram appears to be similar to the B-T graph, because 
as we recall, /i acts as the destabilizer of the boundstate just as T does. The right panel also 
shows the relation /i ~ ±i?^ in the region of small B and the slopes of the plots increase as 
the temperature increases. 

Let us look into the B-n phase diagram a little further. When the temperature is 
high, say, T = 0.3, we have the relation of the critical values of the chemical potential, 
fic{B = 0) < Hc{B = 0.8). In the low temperature regime, say, T = 0.2, we have fidB = 
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0) > fic{B = 0.8). This implies that the fi-T phase diagram, plotted for i? = and 0.8 
must cross somewhere between 0.2 < T < 0.3, that is, the peculiar phenomenon due to the 
fermi sphere must appear as a "cross" in the ^i-T phase diagram. Indeed, the "cross" is 
observed in the fi-T phase diagram shown in Figure [T8l In the low temperature region, we 
observe that the critical chemical potential for the lower B is higher and the higher B in 
this region of temperature is giving smaller critical /x. This counters the intuition that B 
and /X are the stabilizer and the destabilizer of the quark boundstate, respectively. So we 
learn that when the LLL is not dominant due to the fermi sphere, this intuition is false. 

2.6.3 Pure E Background 

In this subsection, we examine the constant electric field background. As remarked in Sec- 
tion [2321 the essential singularities in the effective potential ()64p prevent us from evaluating 
this as it appears. Therefore, we resorted to two approaches, one is the zero temperature 
limit and the other is the weak E expansion. We discuss those in turn. 

Zero Temperature 

The effective potential of this case (without the chemical potential) is given in Equation ()65p 
and we recall that this is always complex for any E ^ 0. The general idea is that we extract 
the dynamical mass from the real part of the effective potential ()67p and check if the 
imaginary part ([71]) is not too large compared to the real part. 

It is straightforward to obtain the tadpole condition (|26p for the real part of the effective 
potential (j67p and numerically solve for the dynamical mass m* for a given E. The result 
is presented in Figure [19] for g = 10 and 20. The diagram indicates that the E field, 
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Figure 19: The plot of the dynamical mass with respect to the E field at T = 0. For the coupling <; = 10 
and 20, we have the second and first order phase transitions, respectively. 

in general, acts as the destabilizer of the quark boundstate, just as the temperature and 
the chemical potential (c./., the pure B background T-m* diagram of Figure [15] and fi-m* 
diagram of Figure [TT]) . When the coupling is not too large (but above the critical coupling 
gc which is finite), the mass continuously goes down to zero, indicating the second order 
phase transition. This second order phase transition was discovered in Reference [21j . 

What has not been noticed is the existence of the first order phase transition at a large 
coupling constant, as shown in the figure for g = 20. The behavior of the effective potential 




29 



is very similar to that of Figure O where the role of the chemical potential in that figure is 
replaced with the E field (and of course T = for the current discussion). For this reason, 
we can expect the very similar result when we include the chemical potential to this zero 
temperature system. 

We now must check if the results drawn from the real part of the effective potential above 
have any physical significance, by examining the imaginary part of the potential. Figure [20l 
shows the on-shell values of the effective potential, that is, it is evaluated at m = m* , with 
respect to the E field. The figure shows that the real part dominates when the E field is 
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Figure 20: The plots of the real and imaginary parts of the difference action at T = and m = m* . 
The curves terminate at the phase transition point. Near the phase transition points, the imaginary part 
strongly dominates over the real part. 



small and away from the critical value. However, when the electric field is near the critical 
value, the imaginary part strongly dominates over the real part. This could indicate that 
the information about the phase transition points is not reliable, due to the breakdown of 
the stable background by the pair creation. 

To obtain further insight into this issue, we recall that Schwinger has worked out the 
rate of the pair creation in Reference [H] and it is roughly given by exp[— vrm^/E']. This 
means that the results from the real part is reliable when > E. Again, this condition is 
not so well satisfied around the point of the phase transitions. Therefore, though there are 
good indications of the first and second order phase transitions, it is not safe to conclude 
so. Most certainly, the values of the critical electric field cannot be trusted. 

As in Reference [21], the imaginary part of the effective potential has not been paid much 
attention, probably due to the NJL model's lack of the confinement. The confinement may 
fix the problem, but then we are dealing with different model and the conclusions from the 
real part still does not seem to correspond to the real world. Nonetheless, we have carried 
out the analysis for both the real and imaginary parts, because we can compare the results 
to the deconfined phase of the SS model. 

Weak E Expansion 

We now discuss the weak E expansion of the T > effective potential (j72p . The expansion 
completely washes out the singularities and the potential is real. One can easily obtain 
the tadpole condition (j26p from the potential (j72p and solve numerically for the dynamical 
mass. Figure [21] shows the result for g = 10, T = 0.3 and up to 0{E^). We observe fairly 
good convergence of the higher terms and the result may be trusted. In particular, the 
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second order phase transition observed at T = appears to persist at finite temperature. 

It is of interest to see if tlie first order transition at tlie zero temperature can also be 
observed at finite temperature. Figure [22] shows the plot of the effective potential at 5 = 20 
and T = 0.3. The 0{E^) effective potential shows the first order phase transition, but the 
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Figure 22: The eflFective potential at T = 0.3 and g = 20. The 0{E^) on the left shows the first order 
phase transition consistent with the T = result, however, the inclusion of the higher order terms [shown 
on the right is up to C'(i5*)] destroys this behavior. 

examination of the higher order terms imply that this may not be trusted. On the right 
panel, the effective potential up to 0{E'^) is shown. We see that the higher order terms are 
not converging and invalidates the result of 0{E'^). Therefore, we cannot conclude if the 
first order phase transition at g = 20 and T = persist in finite temperature. 

When the temperature is lower, the convergence of the higher order terms gets worse. 
In all cases, the potential up to 0{E'^) shows all the expected behavior (as in the previous 
paragraph), but one finds that the higher order terms are behaving completely unphysical 
way. We have expanded the T = effective potential with respect to small E and found 
similar defects in the higher order terms, so the expansion at low temperature is not valid. 
For this reason, we do not attempt to draw the E-T phase diagram. 
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3 Sakai-Sugimoto Model 



We analyze the SS model in this section. We specifically deal with the one-flavor case, but 
the generalization to more than one flavor is straightforward. We flrst discuss the general 
features of the model, including the chiral phase transition and the dynamically generated 
mass of the model. Then, after setting up the model with external parameters, it is further 
examined separately for the conflned and deconflned phases. Finally we present the results 
of the numerical evaluation. 

3.1 Chiral Phase Transition and Dynamical Mass 

The SS model is deflned on the IIA background geometry of D4 branes in the near horizon 
form, 

ds^ = (U/Rf/\ftdXl + {dX,f + hidX^f] + [R/Uf/\r^dU^ + U^d^l) 
$ = $0 + ^ Hu/R) 

C(3) = -2,g~^R^ sin^ Oi sin^ 6*2 cos 6*3^611 A d6'2 A dO^ , (74) 

where dO./^ is the metric for the four-sphere with the coordinates ^j, $ is the dilaton and 
C(3) is the RR potential to which the D4-branes couple magnetically. We also have 

Qs := exp($o) , ■■= TrgsNJl , (75) 

and the "blackening factors" are 

ft = f = l — (Ut/U)^ and /4 = 1 for the high temperature deconflned phase 
/4 = / = 1 — (Ukk/U)'^ and ft = l for the low temperature confined phase . (76) 

The X4 direction is compactifled and at flnite temperature, the Euclidean Xq direction is 
also compactifled as well. To avoid the conical singularity from the compactiflcations, we 
must have the circumference of X4-direction, as 

^, = —{-—] , (77) 



3 \Ukk 

in the low temperature conflned phase and the circumference of Xg-direction, /3 as 



in the high temperature deconflned phase. 

The system with the non-compact X4 direction is of interest as well. In this case, the 
system is always in the deconflned phase and this allows us to discuss the low temperature 
behavior of the model in the deconflned phase. In fact, when we discuss the deconflned 
phase below, we will be cavalier about the conflnement/deconflnement critical temperature. 

Chiral Phase Transition 

The model has the well-known three phases [22] shown in Figure [23l The thin lines in 
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(a) confined %SB (b) deconfined%SB (c) deconfined%S 



Figure 23: Three phases of the SS modeL The coordinate U is the radial coordinate transverse to the D4 
branes. The thin hnes shown in the X4 circles represent the probe D8 branes. The diagram (a) represents 
the low temperature confined phase with broken the chiral symmetry (^SB), (c) is the high temperature 
deconfined phase with the restored chiral symmetry (xS) and (b) is the intermediate deconfined phase with 
XSB. 



the X4 circles are tlie probe D8-branes. Tlie confinement [Diagram (a)] and deconfinement 
[Diagrams (b,c)] phases are determined only through the D4 background, independent of 
the probe D8, and the phase transition occurs at /3 = /34 |22j . 

The U-shape configurations of D8 shown in Diagrams (a,b) represent the chiral symme- 
try breaking in this model. There is a range of temperature and asymptotic separation of 
the D8 branes, denoted by L, so that the background is in the deconfined phase and the 
chiral symmetry is broken, as shown in the diagram (b) of Figure [23l Aharony et al. |22j 
have shown that this intermediate phase exists when L < ^ 0.97(/34/27r) in the decon- 
fined phase. When L > Lc, the intermediate phase does not exist and the chiral phase 
transition and the confinement/deconfinement phase transition occurs simultaneously. The 
chiral phase transition in the deconfined phase is determined by comparing the values of the 
D8 probe action in U- and || -configurations for a given set of parameters. Since the probe 
configuration jumps from one to the other at the phase transition point, this is always a 
first order phase transition. 

String Endpoint Mass (SEP) 

The model in the U-shape configuration exhibits various spectra which are similar to low 
energy QCD, such as mesons. For instance, the fluctuations of the embedding map from 
D8 to the spacetime represent the scalar mesons and the D8 worldvolume gauge fields (of 
non-zero modes in a certain mode expansion) correspond to the vector mesons. Moreover, 
the higher spin mesons are described by the spinning open string of the following configu- 
ration [23]. The endpoints of the string are both attached to the the tip of the U-shaped 
D8-brane, in the different locations of the field theory directions, and they are rotating 
around each other in those directions. The body of the string is draped all the way down 
to the wall and lying there. Figure [241 represents a slice of the string configuration and the 
field theory directions are coming out of the diagram. The "wall" is the tip of the "X4 
cigar" in the confined phase and the horizon in the deconfined phase (see Figure 1231) . 

The energy of the string segment in the U direction (the vertical segment in Figure [M]) 
is given by 

= ,^,u=—j_ VfJ^^U. (79) 

where U = Ukk and Ut for the confined and deconfined phases, respectively. It is argued 
quite generally in Reference [23] that m* represents the dynamically generated mass of the 
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Figure 24: A slice of the string configuration describing tfie high spin mesons. The location U = Uo is 
the tip of the U-shaped D8 and U — Ukk for the confined phase and U = Ut in the deconfined phase. The 
string is also stretching out of the figure into the field theory directions and coming back up again to the 
D8 tip. 



and vector mesons increase monotonically with the energy of the vertical segment and when 
the vertical part is long enough, the masses increase linearly. Also it is pointed out that the 
model (without a tachyon potential) is believed to describe a QCD-like system with zero 
current quark mass, because the pions are massless. This leaves the only candidate, the 
dynamically generated quark mass. This mass, however, is different from the constituent 
mass and is called "string endpoint mass (SEP)"0 The difference mainly comes from the 
energy of the string segment that is horizontal to the wall. For this reason, we refer this 
mass as the dynamical mass or SEP. 

One might argue that the meson spectrum and the dynamical mass would not make 
sense in the deconfined phase. However, we have the fact that the SS model exhibits the 
chiral phase transition in the deconfined phase, and moreover, there are evidences from the 
lattice simulations that indicate the existence of the quark boundstate in QCD, significantly 
above the deconfinement temperature 00]. See Reference [24j for the holographic study of 
mesons in the deconfined phase. 

In what follows, we will be interested in the dynamical mass described here and we 
examine it for both in the confined and deconfined phases. 

3.2 External Parameters 

Constant NSNS Field 

The aim of this section is to investigate the chiral phase transitions and the dynamical 
masses under the influence of external parameters, such as temperature. One of the external 
parameters in which we are interested is the background constant electromagnetic field and 

" The arguments made in Reference is especially reliable for a very high spin mass, because the 
semi-classical treatment of the string is good in this regime. However, the analysis of the Regge trajectory 
carried out in the reference indicates that the description of the mass is fairly good even for the low spin 
mesons. 

Private communication with J. Sonnenschein. 




One evidence shown in the reference is that the masses of the massive pseudo-scalar 
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we would like to introduce this into the SS model. The most natural thing to do is to 
introduce the background gauge field of diagonal U{l)v that comes from U{1) x U{1) of 
D8 and D8 probes. This, in fact, behaves exactly like the electromagnetic field for our 
one-flavor model0 Now, recall that the gauge invariant quantity in the D-brane action is 
not the field strength tensor of the gauge field, but the combination of that with the NSNS 
-B-field (-Bns)) pulled back on the brane with respect to the embedding map of the D-brane. 
Therefore, we can freely trade the background gauge field to i?NS) up to the embedding 
map and the factor of 2Tra' . As we will describe momentarily, we will choose the standard 
simple embedding map of the D8 brane so that the introduction of the constant background 
electromagnetic field is equivalent to the constant i?NS in the IIA supergravity background. 
Notice that the constant Sns does not affect the supergravity background. 

Let us then consider which components of -Bns we should turn on. We are interested 
in the fields that have direct interpretation in the dual field theory. This implies that we 
should have at least one component in the field theory direction. We also avoid turning on 
(5ns )m) because according to the general holographic dictionary, they source the boundary 
currents in Xj-direction and we do not have such counterpart in the NJL model. [We will 
discuss the exceptional case (-Bns)o(7 in a moment.] There are interesting possibilities to 
turn on (i?Ns)Aj4 and (-BNs)4r/- However, these terms complicates the probe action and 
more importantly, the field theory interpretation of those quantities are not clear, so in the 
present work, we will not discuss these possibilities. We are then left with 

Ei := (-BNs)oi and Bi := eijk{B^s)jk , (80) 

where €123 = +1 for the convention and we deal with these background fields in the following. 

Before setting out for the further analysis, we must mention the topological constraints 
imposed on the constant i?NS- As can be seen in Figure [23l the SS model possesses the 
degenerate point in the background geometry, and any constant tensor (of rank> 0) cannot 
make sense as a field at the degenerate point. This severely restricts the possible constant 
Sns and the situation is summarized in Table [H together with the excluded cases for the 
completeness. In particular, notice that the electric fields Ei are not at all allowed in the 
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Table 1: The topological constraints on the constant _BNS-field. The top row is the components of BNS-field. 

deconfined phase. As we will see, the free energy of the system in this case is always imag- 
inary for any nonzero constant electric field and this is precisely due to the violation of the 
topological constraint. The reader who looked at the previous section for the NJL model 
would recall the same imaginary behavior in the effective potential. There, we carried out 
the analysis observing the real and imaginary parts separately. Thus, we do here the same 
for the sake of the comparison and we assume that the imaginary part has the same physical 
interpretation as in NJL, namely, the instability due to the pair creation of the quark pair El 

As remarked in Reference [55], the quarks have different couplings to the electromagnetic field for 
Nf > f , but even in this case, we should expect that the U{l)v gauge field acts similar to the real background 
electromagnetic field. 

We do not induce the current a la Karch and O'Bannon [JT], because we do not have such counterpart 
in the NJL model. See Reference [26] for the treatment with the current. 
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Chemical Potential 

In addition to the constant electromagnetic field, we include the quark chemical potential 
and see how the system respond to this. For this, let us consider the (0, ?7)-component 
of the i?NS field. The AdS/CFT dictionary tells us that this sources the charge density 
of the boundary theory. Recall that a chemical potential can be thought to be the time- 
component of the gauge field that minimally couples to the conjugate charge density. Thus, 
this component of i?NS-field should be directly related to the quark chemical potential of our 
interest. In fact, on the probe D8-brane, this is equivalent to the same component of the field 
strength tensor, indicating that (-Bns)o!7 is the [/-derivative of the chemical potential. The 
treatment of the chemical potential in Reference [25] actually starts by including the time 
component of the world volume gauge field and we have just argued that this is equivalent 
to the introduction of (i?Ns)oc/- Thus, in addition to the components ([80|) . we include 

Eu{U) := B^souiU) , (81) 

where following Reference [25], we assume the dependence of this component on U and this 
function must be determined by solving the equations of motion that follow from the probe 
action. Finally, notice that since 

d[Eu{U)dXoAdU]=0, (82) 
this does not affect the supergravity background. 

3.3 Probe Action and Equations of Motion 

Embedding Map 

Writing down a D-brane action requires the information on how to embed the D-brane 
in the spacetime. We adopt the standard embedding map of the original paper [5]. To 
straighten out the notations and conventions, we explicitly write out the map. 
Let the spacetime coordinates be 

spacetime: X^, /i = , 1 , 2 , 3 , 4 , u , ,02 ,63 , ^4 , (83) 

and let the D8 worldvolume coordinates be 

B8: Xa, a = ,1 ,2 ,3 ,14,6*1 ,6*2 ,6*3 ,6*4 . (84) 

Note that a = 4 is missing. With the slight abuse of notation, let X denote the embedding 
map 

X : D8 spacetime , (85) 
and following the standard practice, this is the identity map, except 

Xi{Xa) = X4,{Xu) , (86) 

where x^i^Xu) is the profile function. 

We define the metric on the D8-brane as the pullback of the spacetime metric two-tensor, 
with respect to the embedding map, 

5D8 := X*g,t , (87) 
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and this yields 



ds 
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u^dn 
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ftdxl + {dxif + ^x'ih + 

where we have renamed and x'^ := dx^/du. Parenthetically, we also have 

ft = f = l — (ut/u)^ and /4 = 1 for the high temperature deconfined phase 
fi = f = l — (tiKK/'ii)^ and ft = ^ for the low temperature confined phase , 

where the spacetime coordinate U in Equations (j76p is changed to the worldvolume coor- 
dinate u. In addition, we pullback i?NS onto the D8 worldvolume, i?NS '■= X*B]<fs and we 
have 
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(90) 



where we have dropped i?2,3 and B2, because we will be interested in the system with 
B {A. OY \\) E and for this purpose, the components above are sufficient. We have also 
extracted the factor of i and redefined the real Ei^u- 

Probe Action 

For our background C(3) and -BnSi we do not have the Chern-Simons coupling of the D8- 
brane, so the probe action consists only of the DBI part 



S-DBi = 2T8 / d^xe J det(5tD8 + -Bns) , 



(91) 



where Tg is the brane tension and the factor of 2 comes from the D8 and D8 contributions. 
We define the "Lagrangian" C through the relation 



where 



det(5D8 + Bss) = Qs \^\\ , 



1^4! := sin^ 9i sin^ 62 sin ^3 



(92) 



(93) 



According to the discussion made in Section [2m the cases B \\ E (B3 = 0) and B _L 
E (Bi = 0) correspond to the U{1)a anomalous and non-anomalous boundary theory, 
respectively. Though we are interested in the non-anomalous case at the end, we carry 
out the analysis for both cases up until the numerical evaluation, because it is easy to 
generalize to Nf > 1 and we can observe some general similarities and dissimilarities to the 
NJL model. Examining separately for those cases, we obtain for the anomalous case 



{{R/uff^' + fx'^^}Bi£i - {R/ufBiEl 



1/2 



(94) 
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and for the non-anomalous case 



r — 4 

''-NAnom — U 



{{R/uff^' + fx',''} A - {R/ufBsEl 



1/2 



(95) 



where we have defined 



Bi:=l + {R/ufBf , £i:=l- {R/uff^'El 

and A:=l + {R/uf{Bl - fi^Ef) . (96) 
We see that Cs depend only on the coordinate u, so the action can be written as 

Sdbi = 2uj,V^Tsg;'P J duC , (97) 



where (3 is the integration of Euclidean time, 



u)i := / ^01^6*2^03^04 1 , and V3 = / d^x 



(98) 



Equations of Motion 

The degrees of freedom in the action appears to be X4{u) and Eu{u). In the D8 worldvolume, 
however, it is more natural to trade i?NS to 27ra'F, that is, we identify (i?Ns)oM = Eu{u) as 
the u derivative of the worldvolume gauge field Aq , and treat the gauge field as the degrees 
of freedom in the action (see Reference [25] )• Thus, we replace as 



Eu{u) ^ /x(u)' := 27raUAo(u) = 2^(^-^ 



lsYd{AoR) 



Rj d{u/R) ' 



(99) 



where we have emphasized that everything is measured in the units of the scale R and is 
dimensionless that involves the factor of {Ig/R)'' in its definition. 

Identifying the degrees of freedom, the equations of motion that follow from the action 

are 



dC , dC 

ci and 



dx'. 



dfi' 



-C2 



(100) 



where ci^2 are the constants of motion and the minus sign on C2 is inserted for the later 
convenience. For /^NAnom of Equation ([95]) . they lead to 



ci =u^fAx', 



{{R/uff^^ + fx','] A - {R/ufBsfi" 



-1/2 



C2 =U%(i?/u)V 



{{R/uff^^ + fx','} A - {R/ufB^ii 



/2 



-1/2 



(101) 



The equations of motion for £Anom of Equation ([M|) can be obtained by the replacements 
A Bi£i and B3 ^ Bi. 
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3.4 Confined Phase 

We now specialize the discussion to the confined phase. We will concentrate on the B J- E 
non-anomalous case, Equation ([95]) . because the anomalous case, Equation can be 

obtained by the replacement just mentioned. For the clarity of expressions, we set B := B^ 
and E := Ei, accordingly, we have ^ = 1 + {R/ufiB"^ - E'^) and B :=! + {R/ufB'^. Also 
we remember that in this phase, = 1 and /4 = / = 1 — (ukk/'u)^- 



The U-Shape Solution and the Free Energy 

The Xi direction is degenerate at u = ukk in this phase (see Figure [23]) and the D8 and D8 
probes must join somewhere. Following the common practice, we assume that the probes 
join smoothly at u = Uq, forming a U shape, and this results in the condition 

a^K^^o) = . (102) 

The equations of motion (jlOip can be solved for x'^ and /u', then we can use the condition 
()102l) to eliminate the constant ci in favor of uq. By resetting c := C2, we obtain 

,2^ iR/uffoiAo/Bo){4Bo + c\uo/Rf} 

P [{R/uriu^fA - uUoAo) + {cyu^)Wf{A/B) - ulUAo/Bo)}] 

= i^/^f^' (103) 

^ {RMHu^fA - uUoAo) + {cyu^)Wf{A/B) - uMAo/Bo)} ' ^ ^ 

where the quantities Xq are meant to be X{uo). 

It is clear in the equation of fi' that this blows up at u = uq leading to a singular Bns- 
To avoid the singularity, we must set 

c = , (104) 
for this U-shape solution and the solution boils down to 

-'.' = rH-X(4r^-A" -d ^'.0. (105) 

The condition (I104p has the drastic consequence that the system is completely independent 
of the chemical potential /i and in particular, the dynamical mass will be insensitive to the 
chemical potential. 

Consulting with Equation (|97p . we define the (dimensionless) free energy of the system 

by 

F := R-^ / duC , (106) 

J UQ 

and we are going to evaluate the free energy for the solution derived above. Plugging the 
solution (jlOSp in the Lagrangian ()95p . we obtain the on-shell expression 

C = u'(-] '^'/- 1/2^ 1/2 ( 1 _ . (107) 



u J \ u^fA 
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Following Reference [22j, we measure the quantities in the units of uq by introducing^ 

y:=u/uo, vr — R/uq, and 2/kk := -UKK/^to • (108) 
Again, following the same reference, we further set z = y~^, yielding 



F = Ivr'" I' dz.-W6;-i/2^i/2(^i _ z^/^^y'^' . (109) 



From the relation ()77|) . we have 



(47r \ ^ 
— j and (34 := P^/R , (110) 



and we can eliminate i/r in favor of i/kk- We then have 

Jo Ac 



-1/2 



F =l{a-'Plr'/'y^T £ dzz-^'/'f-'/Uy^ (l - .8/3: 
/ =1 - yij^z ,A = l + zplyi^{B^ - E^) , 

/o =1 - y| K , A = 1 + /34'yKK(^' - E^) , (111) 

where we set B := a^^^'^B and E := a~^/'^E. We again note that the anomalous case with 
the Lagrangian can be obtained by the replacement A BiSi. 

The free energy obtained above actually is ill-defined, for it diverges at z = 0, and 
currently it is not known how to regulate this. Still, we can draw general conclusions from 
the free energy. Consider the case with B = 0. Then we see that the electric field may not 
exceed the critical value 

< i^)' ■ (m 

This is the curved space version of the upper limit on the electric field in string theory. 
Beyond this value, the free energy becomes imaginary, indicating the instability of the sys- 
tem. Thus, the electric field acts as a sort of "destabilizer" in the system. On the contrary, 
we note that the B field behaves exactly the opposite; it is the "stabilizer" for the B J- E 
non-anomalous case. In particular, those effects exactly cancel when B = E and the sys- 
tem is completely the same as the one without the external electromagnetic fields. For the 
anomalous B \\ E case, those external fields do not interfere each other and they behave 
more or less independently. Later, we will compare these behaviors to those of the NJL 
model. 



Inter-Brane Distance and SEP 

The SS model has another parameter L which is the asymptotic separation of the D8 and 
D8. This is related to the parameter yxK := uyj^/uq and we prefer L to yxK as the 
fundamental parameter of the theory. 

^■^ This scaling by uq is a temporary prescription and as mentioned earlier, we will measure all physical 
quantities in the units of R at the end. 
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For the relation between those parameters, we start with the equation 

r f'oo 

L= dx4 = 2 x^'du , (113) 



J Jua 

where the factor of 2 accounts for D8 and D8. Using the solution (jlOSp and carrying out 
the changes of variables u ^ y ^ z as before, we obtain 

I :=L//34 

l'\,f-l(,-lf± _ ,^I^Y" , (114) 



where we have defined L := L/R and £ has the range < £ < 1/2. 

Now the dynamical mass (SEP) discussed in Section 13.11 can be written as 



= l^i^o/R) C fivY^'^dy . (115) 



We notice that the mass is totally independent of the temperature and the chemical poten- 
tial. It turns out that the y integral can be carried out to yield a closed form as a function 
of yKK- Trading the parameter uq/R to (3^ and t/kk using the relation (jllUp . we get 

m* :=27rf^^ {m*R) 



=a/34-2{ - V^r(2/3)/r(l/6) +yKK2^i(-lAl/2,2/3,yiK)} • (116) 

This dimensionless mass m* together with Equation (jll4p is ready for the numerical eval- 
uation. 



3.5 Deconfined Phase 

We now proceed to examine the deconfined phase. Strictly speaking, the discussion here 
would make sense only above the confinement /deconfinement critical temperature. How- 
ever as mentioned earlier, we can always consider a very large /54 or even the non-compact 
X4 model, so we will freely discuss the low temperature behavior of this phase without 
worrying about the critical temperature. Like in the confined case, we will concentrate on 
the B J- E non-anomalous case with the Lagrangian (195p . For the clarity of expressions, 
we set B := B3 and E := Ei, accordingly we have ^ = 1 + {R/u)^{B'^ — f'^E"^) and 
B := 1 + {R/u)^ B"^ . Also we remember that in this phase, /4 = 1 and ft = f = 1 — (ut/u)^. 

The U- and \\-Shape Solutions 

In the deconfined phase, the xq direction is degenerate at u = ut but the X4 direction is not 
(see Figure [23l) . Therefore, the D8 and D8 probes can either join smoothly in the U shape 
or go straight down to the horizon, forming the || shape. Those configurations correspond 
to the broken and unbroken chiral symmetry phases, respectively. 

The derivation of the U-shape solution from the equations of motion (jlUip is nearly 
identical to the confined case and we have 

<' = r'(-Y(4r^-A" ^'.0. (117) 
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Again, the system in the U shape is completely independent of the chemical potential fi 
and in particular, the dynamical mass will be insensitive to the chemical potential. 

For the ||-shape configuration, we have x'^{u) = 0. Thus, the equations of motion (llOip 
yield the || solution 

" ° ^' = B{u^{R/u)^B + c-} ' ^''^^ 

where we set c := C2- The anomalous B \\ E solution can be obtained by the replacements 
A Bi£i and B Bi. 

We notice that for any finite E, there is a factor f~^E'^ inside A (or £) and this diverges 
at ti = UT- To avoid this, we must conclude that we need c = 0, i.e., //' = for any 
E ^ 0. Moreover, since fi must vanish at the horizon, we get /i = 0. This strong conclusion 
has a clear physical interpretation. In the deconfined and chirally symmetric phase, quarks 
and anti-quarks are not bound together and the chemical potential causes the unbalance 
in the charge density. Thus, application of the external electric field immediately triggers 
the instability. Note that the key factor in this instability is f~^E'^, which is the direct 
consequence from the violation of the topological constraint discussed in Section 13. 2[ 

The Difference Action 

Given the two solutions, we must determine the energetically preferred configuration for a 
set of external parameters. For this purpose, we compute the difference action with respect 
to the solutions. As before, we define the dimensionless free energy 

F ■= R-^ J duC , (119) 

where C is given in Equation (j95p . Evaluating the Lagrangian for the solutions (|117p and 
(fTTHI) . we get 



Cu =u 



4 / R\ 



u 



3/2 / c'^ \ 



uj V u^iR/u)^B 

Now the difference of the free energies is 



(120) 



oo 



AF = R-'' / duCu - R / du£^^ , (121) 

so that when this is negative (resp. positive), the U shape (resp. || shape) is preferred. For 
the numerical handiness, we rewrite this as 



u^{R/ufB 



Lr r\r) V^u^iR/ufB) ■ ^^''^ 
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We carry out the sequence of the variable changes, u ^ y ^ z, as was done in Section [32 
Assuming the relation 

— j and (3:=[5/R, 
that follows from Equation (|78p. we have the expression of AF in terms of ht as 



(123) 



,8/3 



dzz-^^/^A^lH 1 



+ z 



5/3 



c 



-1/2 



-1/2 



l + Z 



5/3 



-1/2 



y^z 



fo = i-yT 



where C := c^/uq and we have 

/ = ! 

and setting B := a'^^'^B and E := a^'^^'^E, 

A=l + zPS^B^ - f-^E^) , B = l + zPS^B^ , £ = 1- zP^y^f-^E^ 
Ao =1 + P'^vUB^ - fo'E^) ,Bo = l+ P'v^tB^ ,£o = 1- PWrfo'^^ ■ 



(124) 



(125) 



(126) 



As noted before, because of the factor f^^E"^ in A and £, the integrand of the second line 
in Equation ()124p is divergent at z = y^^, however, it turns out that the integral is finite. 
When we consider the pure E background, we find that the difference free energy is complex 
and the imaginary part comes only from the || configuration. But since we cannot define 
the action without taking the difference, we treat the complex value as representing the 
nature of the whole system. 

Chemical Potential, Inter-Brane Distance and SEP 

The II configuration without the E field has the non-trivial /i' as shown in Equation pisp . 
We would like to set the asymptotic value of fi to p,, which should correspond to the chemical 
potential of the dual field theory. We must also require that fJ-{uT) = 0. Therefore, 



1^ 



du 



Ac^ 



13{u^{R/u)^13 + c^} 



1/2 



dzz-'/^ 



AC 



e{a-3/36y3^ + C7z5/3} 



1/2 



(127) 



The asymptotic brane distance for the ||-shape configuration is a parameter on which 
no other parameters depend. That of the U-shape configuration can be worked out in the 
similar manner as in the confined case and we get 



L :=L/R 

_2 
~3 



dzf 



-1/2 



1 fA 



foAo 



,5/3 



-1/2 



(128) 
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The dynamical mass is simply proportional to the distance {uq — ut), because we have 
V9009UU = 1, and we get 



m* :=27rj^-^J {m*R) 

=ap-\yT' - 1) . (129) 

The expressions for AF, fl, L and m* obtained in this subsection are ready for the nu- 
merical evaluation and they suffice to determine the dynamical mass and the phase structure 
of the deconfined phase. 

Zero Temperature 

For the sake of the comparison with the NJL model, it is interesting and important to 
consider the deconfined phase in the low temperature regime, including at zero temperature. 
At zero temperature, the parameters of the theory change slightly. We must replace the 
"blackening factors" as /, /o ^ 1 and we have the range < n < oo where ut is no longer 
the parameter of the theory. We then have 

^ := 1 + zy^a^iB"^ - E'^) , B := 1 + zy^a^B"^ ,£:=!- zy\a^E'^ , (130) 

where we have extracted the factor of a by B^ = a^^a^^B"^) = a^B^ and similarly for E, to 
be consistent with the finite temperature parameters. 

The difference free energy shown in Equation (I124|) does not change, except that the 
parameter yn is the parameter of the theory itself and not expressed in terms of yT, and 
the replacement (JT30]) is assumed. Examination of the difference free energy shows that the 
U-shape configuration is preferred to the || configuration at zero temperature. 

The chemical potential and the inter-brane distance for T = take the forms 



- 1 -1 



1/2 

dzz-^l^ 







AC 



B{y%B + Cz^/^] 



It is also easy to show that the mass becomes 



r^*=2^(^|) {m*R)=y-^\ (132) 
3.6 Numerical Evaluation 

We now present the numerical results, separately for the confined and deconfined phases. 
In this section, we assume all the quantities are measured in the units of R and we are going 
to omit the hats {X) on the parameters. 

3.6.1 Confined Phase 

We have seen in Section 13.41 that the system is independent of the temperature and the 
chemical potential, except that this phase is defined below the critical temperature of the 
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confinement/deconfinement phase transition. Therefore, the only parameters are the exter- 
nal electromagnetic field and the inter-brane distance. As for the numerical evaluation, the 
relevant equations are £ := L / of (I114p and m* of ()116p . 

We first examine the properties related to the parameter i without the background 
electromagnetic field. From Equation (I114p . we can plot i with respect to yxK and it is 
shown in Figure [25l Recall that the ranges of i and yxK are (0, 1/2) and (0, 1), respectively, 

B=0=£, /?4=1 




0.0 0.2 0.4 0.6 0.8 1.0 

yKK 



Figure 25: The plot of I with respect to the variable j/kk without electromagnetic background and (34,-1. 
This shows that the location of D8 and D8 can be anywhere in the X4 circle. 

and the figure shows that the parameter i takes all the possible values. We remark that the 
lower bound 2 of is the geometrical constraint and not the bound determined by the 
dynamics (like the NJL coupling g^)- We will further discuss the comparison between 
and g later in Section HI 

Figure [26] shows the response of SEP with respect to £^^. The left panel shows the 



B=0=£, /J4=l B=0=E. /34=1 




50 100 150 2.00 2.02 2.04 2.06 2.08 2.10 2.12 2.14 

r' 



Figure 26: The plots of mass m* against (2 < < cxa). The left panel implies the relation 
m* ~ (^~^)'^ for large values of and the right panel shows that the relation m* ~ is fairly well for 
small values of £~ . 

relation m* ~ (^~^)^ for the range of large away from the lower bound and the right 
panel shows that near the lower bound 2, we have the linear relation m* ~ £~^. 
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Pure B Background 

When the B field is turned on, the vkk-^ graph appears essentially the same as Figure [25] 
and the parameter (. takes all the possible value < ^ < 1/2. The response of the mass 
with respect to B is shown in Figure [271 The graph shows the general tendency that B acts 




Figure 27: The SEP against the B field. We have defined mo := m*{B = 0) and notice that the absolute 
value m* is much larger for £ = 0.1. The right panel shows the relation m* ~ for small values of B 

as the stabilizer of the boundstate. The right diagram reveals the relation m* ~ B^ in the 
region of small B. 

Pure E Background 

When we turn on the E field, the ukk-^ relation becomes more complicated. The left panel 




0.8 1.0 




f = 0.15. ;ii„= 19.85 \' 
r = 0.25, »!„ = 7.976 
r = 0.35, nil, = 3.856 



0.4 0.6 0.8 1.0 



Figure 28: The left panel shows that for E < 1, £ takes aU allowed values < ^ < 1/2, and for E > 1, £ 
has upper limit, implying that the D8 and D8 may not be separated too far. Notice in this case that there 
are two solutions of j/kk for a given £. The critical value E = 1 shown in this figure depends on the value 
of /34, and in general, the critical value is given at [fS^E) = 1. At this value, the curve approaches the point 
{yKK,£) = (1, 1/4). The right panel shows the response of the dynamical masses with respect to E. The 
parameter mo is the value of m* at J? = 0. For 1/4 < £ < 1/2, the dynamical mass goes smoothly to zero 
at £■ = 1 and beyond which the system becomes unstable. For £ < 1/4, the system becomes unstable before 
the mass becomes zero. 

of Figure [25] shows that there is a critical value {(3fE) = 1, above which i takes all possible 
values and below, there is an upper limit which is lower than 1/2. Therefore, for sufficiently 
large value of (PfE), the D8 probe system becomes ill-defined. Since the supergravity 
background is in the confined phase, the physical implication of this instability is unclear. 



46 



Notice that when the E field is below the critical value, there can be two solutions of 
yKK for a given I. One can show along the analysis of Reference [5] that the lower ?/kk 
solution is the stable configuration and the other is unstable. Hence, we always pick the 
lower yxK solutions to obtain the dynamical masses. The right panel of Figure [28] shows 
how the dynamical mass respond to the E field. As expected, the graph shows the general 
tendency of E as the destabilizer. As discussed in the previous paragraph, there is a critical 
value of E and the curves in the diagram terminate there. 

3.6.2 Deconfined Phase 

As mentioned in Section 13.51 we will be cavalier about the confinement /deconfinement 
phase transition temperature and discuss this phase in all temperature range, including 
zero temperature. 

Let us first examine the range of parameter L without the electromagnetic background. 
Unlike the confined phase, L in the deconfined phase depends on the temperature, as shown 
in Equation (|128p and the zero temperature case is given in Equation (|131|) . Figure [JS shows 
the plots of L for different temperatures. The left panel of the zero temperature plot shows 

T=0, B=0=E B=a=E 
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Figure 29; The plots of the inter-brane distance L in the U-shape configuration. The left panel is for 
T = and this case is assumed to be the non-compact model. There is no upper limit on L and it increases 
as y^yij. The right panel is plotted for finite temperatures. Any finite temperature induces the upper limit 
on L, i.e., D8-D8 inter-distance may not be too large to form the U shape. 

that there is no upper limit on L and it increases as ^/yn- This is a manifestation of the fact 
that at zero temperature, the U-shape configuration is always preferred to the parallel one, 
as mentioned in Section 13. 5i The right panel shows that finite temperature introduces the 
upper limit on L. For a fixed value of L, a sufficiently high temperature yields no solution 
for ut (e.g., at L = 1.0 and T = 0.2 of the diagram), implying that the system is actually 
in the chirally symmetric ||-shape configuration. Here again, the temperature plays the 
role of the destabilizer of the quark boundstate. Meanwhile, a sufficiently low temperature 
yields two solutions for yx {e-g-, at L = 1.0 and T = 0.1 of the diagram). As discussed in 
Section 13.6.11 the smaller solution is the stable configuration and we will always pick this 
one in what follows. 

Before introducing the background electromagnetic field, we would like to examine the 
response of SEP with respect to L~^. One can obtain the mass as a function of from 
Equations ()131|132p for T = 0, and for finite temperature, from Equations ()128|129p . The 
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Figure 30: The plots of m* against For T = (upper left corner), the relation m* ~ holds 

for all values of . At finite temperature, the relation is more complicated as shown in the upper right 
corner. The bottom two indicate the relation m* ~ {L~^y^^ in the region of small L~^, but slightly away 
from the critical values. 



results are shown in Figure [30l The upper left panel of the figure shows the T = case and 
one can show analytically for this case that m* increases quadratically with respect to L~^. 
The upper right panel shows the finite temperature cases and these are more complicated 
than the zero temperature. We have observed in Figure [29] that there is an upper bound 
on L for a given temperature. This translates to the lower bound on and the plots 
of the upper right panel terminates at those critical values of L^^. The diagrams of the 
bottom row are plotted for the region near the critical values of and show the fairly 
good relation m* ~ V slightly away from the critical values. 

Pure B Background 

We now consider the effects of the constant background B field. At T = without the B 
field, we have seen that there is no upper limit on L, and finite i?-field does not change 
this conclusion. Neither does it change the general observation that finite temperature 
introduces an upper limit on L. In short, L behaves similar to Figure [29l except that the 
larger B allows the larger separation L for a given ur^t- This, of course, implies that the 
B field acts as the stabilizer of the quark boundstate. 

In the deconfined phase, the system exhibits the chiral phase transition from the U- 
shape to the ||-shape configuration, which cannot take place in the confined phase. Unlike 
the way we presented the NJL results, we first show the phase diagram before the response 
of the dynamical mass with respect to the external parameters. This is because we must 
consult with the chiral phase diagrams to determine exactly where the mass drops to zero 
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Fi gurc 31; The B-T phase diagram at L = 0.1. The right panel shows the region of small B and the 
parameter To is the critical temperature at B = 0. It shows the relation T ~ B^ in the region of small B. 
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Figure 32: The B-^ phase diagram at L = 0.1. The right panel shows the region of small B and the 
parameter /io is the critical value of the chemical potential at B = 0. It shows the relation /i ~ B^ in the 
region of small B. 
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Figure 33: The n-T phase diagram at L = 0.1. The case with B = reproduces the result of Reference 
(the reference plots TL against fiL^). 
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at the phase transition point. 

The phase diagrams can be all worked out from Equations (jl24ll27ll28p and Figures [311 
[32] and [33] show the B-T, B-^. and /x-T phase diagrams, respectively. The transitions are all 
first order as the configuration "jumps" from U- to || -shape. We observe the general trend 
that B acts as the stabilizer and T, ji as the destabilizer. This fact appears as the general 
similarity in the B-T and B-^ graphs. 

Let us now investigate the response of the dynamical mass with respect to the external 
parameters. Figure [M] shows the B-m* graph. As expected, the mass increases as B is 

L=OJ 
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Figure 34: The SEP plotted against B-field. The right panel exhibits the small B region and the parameter 
mo is defined as the mass at _B = 0. The diagram shows the relation m* ~ in the region of small B. 
Typically at high temperature, the lower B region of these curves may sharply drop to zero as the chiral 
phase transition sets in and the point of the drop depends on the temperature and the chemical potential. 

L=OJ 
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Figure 35: The plot of mass against temperature. One must be aware that the low temperature region is 
subject to the confinement/deconfinement phase transition, unless Pa is large enough, and the high temper- 
ature region typically drops to the zero mass sharply due to the chiral phase transition. The point of the 
chiral phase transition depends on the chemical potential as well. 

raised and higher temperature results in the smaller mass. It appears that the dependence 
on the chemical potential is missing, but one must be aware that the chiral phase transition 
can occur when B is small. For instance at T = 1 and /i = 20, the phase diagrams tell 
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us that there is no chkal phase transition for ah values of i?, so the dashed curve of the 
diagram is not modified. However at T = 1 and /U = 40, the chiral phase transition sets in 
at S ~ 4 and the dashed curve for T = 1 in the diagram must be truncated below this value 
and m* sharply drops down to zero at that point. Therefore, the B-m* diagram depends 
on /i implicitly and ^ determines where the mass drops to zero. 

Figure [35] shows the T-m* diagram. The lower temperature region of the diagram is 
subject to the truncation by the confinement /deconfinement phase transition at /3 = /?4. 
Thus, for example, if we have /34 = 1, then the diagram must end at T = 1. The higher 
temperature region can end with the chiral phase transition, by dropping sharply to the 
zero mass. The chiral phase diagrams tell us, for example, that the B = 10 curve must be 
truncated at T = 1 for ~ 57.6 and at T = 1.5 for fj, ~ 32.5. So again, the truncation 
point depends on the value of the chemical potential. 

Pure E Background 

We now discuss the constant E background. As mentioned in Section 13.51 the introduction 
of the chemical potential is not allowed for this case. Moreover, the difference free energy 
A.F in Equation (jl24p is complex and we follow what we have done for the NJL model: we 
treat the real and imaginary parts of the free energy separately. We extract the information 
about the phase diagram from the real part, and we make crude estimate on whether such 
information can be reliable, by checking if the imaginary part is subdominant. We are 
mainly interested in the zero temperature limit, since in the NJL counterpart, we had most 
information at zero temperature. 

From Equations (jl31ll32p for T = 0, we can obtain the dynamical mass as a function 
of E^ and the result is shown in Figure [36l We see that E generally is a destabilizer of 

T=0 

LOO 
0.95 
f 0.90 
0.85 
0.80 
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E 

Figure 36: The SEP as a function of the E field at T = 0. The parameter mo is the m' aX E ^ 0. The 
curves at higher value of E is terminated at the chiral phase transition points. 

the boundstate, as expected, and also the curves always terminate at the first order chiral 
phase transitions. Figure [37] shows the corresponding plots of the free energy in real and 
imaginary parts. This figure shows that the imaginary part is always comparable to the 
real part, especially near the phase transition points. Hence there is a doubt on the validity 
of the masses and the phase transitions shown in Figure [36] However, notice that m*^ 
is considerably larger than E. Recall that in the NJL counterpart, we also estimated the 
stability of the system by considering Schwinger's pair creation rate formula. The rate is 
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Figure 37: The free energy plotted separately for the real and the imaginary parts. The quantity F is the 



right-hand side of Equation (|124p where the factor of {3y1[^) is omitted from the free energy. 



exponentially suppressed by the factor m*'^ /E and this indicates that the system of the SS 
model may have very slow pair creation rate. Therefore, it is reasonable to assume that the 
instability takes place slowly and the information about the mass and the phase transition 
can be meaningful. 

The system at a finite temperature behaves as expected and it acts as a destabilizer. The 
masses at a certain value of E is smaller at finite temperature than the zero temperature 
case and the first order chiral phase transition sets in at a smaller value of E. The E-T phase 
diagram extracted from the real part of the difference free energy is plotted in Figure [ 




Fi gurc 38: The E-T chiral phase diagram obtained from the real part of the free energy. 



4 Comparison 

We now have the data for the comparison between the NJL and SS models. In this section, 
we summarize the results by putting the NJL and SS data side by side. 



4.1 General Observation 

We first make general observations which are not specialized to purely B oi E backgrounds. 
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Let us compare the NJL effective potential and the probe action of the SS model 
in the confined phase with the Lagrangian (jl07p . Both expressions are specific for the 
non-anomalous B J- E background configuration. The former is the expression for T = 
and the latter is independent of the temperature, as long as the system is in the confined 
phase. We have already remarked under Equation (I54p that for the cases with \B\ > \E\ and 
\B\ < \E\, the NJL system behaves as if there is only B- and i?-field, respectively. The case 
with \B\ = \E\ is special and the system is equivalent to the one without the electromagnetic 
background. Also the pure -E-field case can be obtained from the pure S-field case just by 
sending B — > iE, and vice versa. As we can see clearly in Equation ()107p that the same is 
true for the SS model probe action. A short inspection reveals that the zero temperature 
limit of the SS deconfined phase discussed near the end of Section 13.51 also shows the same 
behavior. (One must recall that fi = for E ^ 0.) 

The finite temperature NJL effective potential (j50p shows that the lack of Lorenz invari- 
ance destroys the features mentioned above. Also the probe Lagrangians of the SS model 
in the deconfined phase ()120p lack the properties of the previous paragraph, because of the 
general appearance of the factor f~'^E'^ in A. 

Therefore, we see that even though the actions of the theories are very different, the 
effects of Lorentz invariance (at T = 0) and non-invariance (at T ^ 0) appear the same 
for both theories in a very general sense. However, we remark that for the anomalous case 
B \\ E, apparent similarity is not observed. 

We now look into the relation between the NJL coupling g and the inter-brane distance 
L of the SS model. The possible ranges of those parameters are shown in Figure [3] for the 
NJL coupling (plotted for Tr^^'^/g'^), Figure [25] for £ := L/f3i of the SS model in the confined 
phase and Figure [29] for L := L/R oi the SS model in the deconfined phase. In the absence 
of S-field, the NJL model always has the lower limit on the coupling and below this critical 
coupling, the system is in the chirally symmetric phase. For example, the critical coupling 
for the zero temperature and zero chemical potential is 2tt. The SS model in the confined 
phase also has a lower limit on namely, = 2. However, this lower limit is different 
in nature: it is the geometrical constraint rather than dynamical, and furthermore, at the 
critical point the chiral symmetry is still broken (the configuration is in the U shape) and 
the region < 2 is not defined, of course. The SS model in the deconfined phase at zero 
temperature does not have a lower limit on and the chiral symmetry is always broken. 
However at finite temperature, there is a lower limit on and below this critical value, 
the system is in the symmetric phase. So it seems that the NJL coupling g and of the SS 
model in the confined phase behave rather differently, and g and of the SS deconfined 
phase have general similarity at finite temperature. 

To investigate further into the relation among those parameters, we have shown the 
response of the dynamical mass with respect to the parameters. Figure [6] is the plot for the 
NJL model and Figure [26] is the SS model in the confined phase and Figure [30] is the SS 
model in the deconfined phase. In general, the mass respond very differently, except that 
the higher g or gives the higher mass. However, examination into the region near the 
critical points of these parameters shows that there is a good similarity between g and 
of the SS model in the deconfined phase at finite temperature. In those cases, the mass 
grows as square-root of the parameters. 

The conclusion out of this comparison is that the parameter g of the NJL model and 
of the SS model in deconfined phase at finite temperature correspond to each other fairly 
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well, especially near the critical values, and in all cases, the increase in the parameters result 
in the higher masses. However we observed that they are not in one-to-one correspondence 
and especially at zero temperature, they show the qualitative differences. 

In what follows, we focus on the comparison between the NJL model and the SS model 
in the deconfined phase. This is because we have already seen that the correspondence is 
rather poor for the SS confined phase. Moreover, the confined phase does not have the 
dependence on the temperature nor on the chemical potential, and it does not even have 
the chiral phase transition. These facts make the confined phase rather meaningless in the 
comparison with the NJL model. 

4.2 Pure B Background 

When we turn on the B field at T = 0, the lower limit of the NJL coupling becomes zero 
and any finite temperature brings this to a finite value. (See Figured]). This makes the 
correspondence of the parameter g of NJL and of SS even better. 

Let us proceed to consider the dynamical mass in the presence of the background B 
field. Compare the B-m* diagram of the NJL model at /i = in Figure [12] and that of the 
deconfined SS model in Figure [M] We see the striking similarity in how the masses respond 
to the B field between the models. In the small B region, they both show the relation 
m* ~ B^ and this fact is more than just the general tendency of B as the stabilizer of the 
boundstate. The figures also show that the ways the masses in the models respond to the 
temperature are also similar. 

We have, however, only looked at the /x = case so far. The finite chemical potential 
reveals the glaring differences between the models. The NJL model exhibits a very com- 
plicated behavior at finite chemical potential due to the de Haas-van Alphen oscillation 
and the destabilizing effect of B at large /i, as shown in Figure [13] (see also Figure [T^ . 
In contrast, the Figure [34] of the SS model does not change the shape at finite chemical 
potential, except that the lower B part of the curves can be truncated to zero mass as the 
chiral phase transition sets in. The dynamical mass of the SS model, SEP, is independent 
of the chemical potential (except at the chiral phase transition point), and this is why the 
li-m* diagram of Figure [11] for the NJL model does not have the counterpart in the SS 
model. 

We can also see the similarities and dissimilarities due to the chemical potential in the 
phase diagrams. For example, we can look at the B-T phase diagrams of Figure [TBI for the 
NJL model and Figure [31] for the SS model. When the value of the chemical potential is 
small, the NJL diagram appears very similar to the SS diagram. However, the "dent" in 
Figure [16] of NJL is due to the destabilization effect of B at large /i and this is not observed 
in the NJL diagram in Figure [311 Accordingly, the "cross" of the T-/x phase diagram in 
Figure [18] does not have the counterpart in the SS diagram in Figure [33] 

Therefore, the models in the background B field behave very similar to each other when 
the chemical potential is zero or small, but as soon as the chemical potential becomes 
large enough to trigger the oscillation and the B-&e\d destabilization, the models behave 
differently. 

We also point out that the order of the chiral phase transition in the SS model is, by 
definition, all first order, because it is defined to be the point where the configuration jumps 
from U to II or the other way around. Meanwhile in the NJL model, the transition is in 
general second order, as long as the chemical potential is low enough. When the chemical 
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potential becomes large, there is a mixture of the first and second order phase transitions 
in the phase diagrams. 

4.3 Pure E Background 

Turning on the background E field introduces the lower bound for both g and even at 
T = 0, and the effect is similar to the finite temperature discussed in Section 14. 1[ Thus in 
this case, the qualitative agreement between g and is good at all temperature. 

Just like the temperature, the background electric field acts as the destabilizer in the 
both theories and this can be seen in E-ni* graph of Figure [19] for the NJL model and 
Figure [36] for the SS model. The NJL model shows both first and second order phase 
transitions depending on the coupling, and the SS model, by definition, shows only the first 
order transition. We noted that the m*^ to E ratio is very large for the SS model and this 
implies the slow rate of the pair creation. Therefore, the phase diagram of the SS model is 
more reliable than that of the NJL model. 

In both models, the free energy is complex for any finite E, implying that the models 
are unstable in the E field background. We have seen that the real and imaginary parts 
are comparable to each other in both models. They are shown in Figure [20] for NJL and 
Figure [37] for SS. 

5 Conclusion 

We have seen good agreements between the NJL model and the deconfined phase of the 
SS model. Because the NJL model is nonrenormalizable, the cutoff must be introduced 
by hand and this makes the energy scale of the model unclear. One popularly advocated 
interpretation is that the model describes the very low energy regime of QCD and the cutoff 
is around Aqcd- Another interpretation is that since the model lacks the confinement, it 
describes QCD in the deconfined phase, but not too far from the phase transition point. As 
far as the comparison with the SS model is concerned, our study strongly implies the latter 
interpretation. 

We have also seen the general discrepancy of the phase transition orders, and the most 
prominent differences between the models are the de Haas- van Alphen effect and the desta- 
bilizing effect of B at large //, observed in NJL but not in SS. The cause of the differences 
can be reduced to the general fact that the SS model in the U configuration is totally inde- 
pendent of the chemical potential. Physically, if a model exhibits hadrons as boundstates 
of the quarks, the spectrum should certainly depend on the chemical potential. The SS 
model of the form we studied in this work does not show this behavior. Moreover, the lack 
of the de Haas- van Alphen and the destabilizing B effects imply that the SS model fails to 
incorporate the fermi sphere. Therefore, it appears that we need extra ingredients to the 
SS model, such as different brane configurations (which corresponds to additional phases) 
and extra terms in the probe action. Further investigation in this direction would deepen 
the understanding of the model's universality. 

In this work, we have concentrated on the case with Nf = 1 and compared the dynamical 
mass of NJL and SEP of SS. It is of great interest to generalize the comparison to Nf > 1 
and investigate the hadron spectrum, even quantitatively, and the effects of other features, 
such as the U{1)a_ breaking terms in both models. 
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A Path Integrals of the Worldline Fields 

We are going to evaluate the path integrals Ib and Ip defined in the beginning of Subsec- 
tion 12. 5i The derivation is not original and can be found elsewhere. We largely collect and 
reproduce the relevant parts of Reference [11] for the completeness, and specialize to the 
finite temperature system at the end. 

We first consider the worldline fermionic part Ip- We split the field ip into the part that 
satisfies the equation of motion ^ (classical part) and the fluctuation part r] as 

V'^(t) =M/^(r)+r?^(r) , (133) 
where they both satisfy the anti-periodic condition. The classical part ^(t) satisfy 

{dr - 2iF)^,^4T) = . (134) 
We can easily solve the equations of motion to get 

*^(r) = {e''^n,C. > (135) 
where are some constants. Since we have 

^^{0) = C^ and ^^{s) = {e^^^')^^C,, (136) 
the anti-periodic condition ^'(0) = — ^(s) implies that 

(l + e2^^«)^^C7, = 0. (137) 
Assuming that the first factor is invertible, we get = 0. Hence, we conclude tha10 

^/.(r) = . (138) 

We are therefore left with 

If = ^PV^(r) exp[-^ J" drr] ■ [Or - 2iF) ■ r/] = (det^[9, - 2iF]fl^ 

= (detA[a.])^/'(detA[l - 2iFd;^]fl^ , (139) 

where the subscript of the determinant reminds us that the field is anti-periodic on the 
time circle and this accordingly affects the momentum. Notice that since the field r/(r) is 



As in Reference P3], the action must be augmented by a surface term in general. However, in our 
case, the classical part vanishes and the surface term does not affect our discussion, hence omitted from the 
beginning. 
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real, the functional determinants have the power of 1/2. Following Strassler [8j, we adopt 
the normalization 

4 

{detA[dr])'/^ = J2{a\a)=4. (140) 

Q=l 

For the other factor, we quote the result from Reference |14j . 

(detA[l - 2iFd-^])^/^ = (detL[cos(Fs)])i/2 ^ (^4^) 

where det^, is the determinant with respect to the Lorentz structure. Putting those results 
together, we obtain the expression in Equation (j47p . 

Let us now deal with the bosonic part Ib- As it always happens in a finite temperature 
field theory, the zero winding mode of bosonic sector needs a special care. As shown in 
Figure [21 this mode has closed paths. Among the paths, there are constant functions 
x{t) = xq and this makes the functional determinant ill-defined. Therefore, the constants 
xq should be separated from the path integral and they all should be summed (or rather 
integrated) over. However, this integral is precisely the one extracted in the definition of 
the effective potential in Equation ()24p . So in practice, we do not have to worry about these 
constant paths and we can proceed by assuming that they are absent. 

By virtue of the Schwinger-Fock gauge, the Lagrangian CxSF is quadratic and we rewrite 
it as 

rx{s)=y ^ f-s 

Ib =M / Vx{t) exp[-- / drx -{-8^ + 2ieFdr) ■ x 

Jx{o)=x 4 Jq 

- ^{x{s) ■ x{s) - x{0) ■ x{0)}] , (142) 

where we kept the surface terms, because they will have the contribution from the classical 
paths. We split the field into the classical part X{t) and the fluctuation z{t) as 

x(r) = Xir) + z{t) , (143) 

where they satisfy 

(-92 + 2ieFdr)^,uX,{T) = , (144) 

and 



X{0) = X , X{s) = y and z{0) = = z{s) . (145) 
The equations of motion ()144p with the boundary conditions can be easily solved and yield 

X(r) = x + i^^^4^(^-^)- (146) 

^ ^ 1 - exp(2ieFs) ^ ^ ^ 

This gives the contribution to the surface terms of Ib, 

1 1 

--{x{s) • x{s) - x(0) • i(0)} = --(y - x) • eFcot(eFs) ■ {y - x) . (147) 
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The path integral now is over the fluctuation z{t) and this gives 

Af{det'p[-d^ + 2ieFdr])-^/^ = Af{det'p[-d^]'^/^){det'p[l - 2ieFd;Y^/^) 
sin(e-Fs) 



eFs 



-1/2 



(148) 



where the prime on the determinant indicates that we omit the constant paths, we have 
used Equation (j35p and the evaluation of the last determinant in the first line is quoted 
from Reference [14j. We thus have 



Ib = (47rs) 



-d/2 



exp 



-{y — x) ■ eF cot (eFs) ■ {y — x) 



deti 



sin(eFs) 
eFs 



-1/2 



Now according to the finite temperature prescription (|37|) . we should set 

y = X + n[3xQ , 

where xq is the unit vector in the time direction. We then have 



Ip = (47rs)-^/2exp [ - (n/3/2)2{eF cot(eFs)}oo] ( deti 



sin(eFs) 
eFs 



^1/2 



(149) 



(150) 



(151) 



B Computation of the Functions of F, 



1/2 

In this appendix, we compute {det^ [(Fs) cot(Fs)] } and {(Fs) cot(Fs)}oo- The former 
can be evaluated along the method of Schwinger [15] and the latter is difficult and we adopt 
the method of Reference [38] (see also [39]). Both discussions are recast to the Euclidean 
version in below. 

We have the Euclidean version of the field strength 
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(152) 



and in the four spacetime dimension, we have the dual of the same rank 
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(153) 



where we have used the convention £0123 = +1- We then have the invariants 



--F F 



2^ 



F^) and Q 



' ^Efj^uFfj^u 



E ■ B 



(154) 
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and 



We first calculate the determinant. The important identities are 



(155) 



(156) 



Using those, we derive the eigenvalues of F^j^iy in terms of and Q. Let tp^ and / be the 
eigenvector and eigenvalue of Ffj_^, respectively. We then have 

Fij,uFuXtpx = fF^ytpy = iGipf, , 

and this implies that 

F^.u'fpu = if~^Gfpfi ■ 

By iteration, we have 

Ff,uF,xi^x = fijf, and F^.F.aV'a = -/"'^Vm ■ 
Then the identity (jl55p gives the eigenvalue equation 

f + 2Tf -g = o, 

with the solutions ±/-|- (signs are arbitrary, i.e., 4 values) where 



(157) 
(158) 
(159) 
(160) 

(161) 



Note that any function of a matrix, $(-F), is defined as the formal power series. There- 
fore, the argument of a determinant can be diagonalized and the determinant reduces to 
the product of the functions Ili^{fi), where fi are the eigenvalues of F. We thus have 



detL[{Fs)cotiFs)] = {s2/+/-Cot(/+s)cot(/_s)}' . 

Since we have 

iuf-f = -g' , 

we immediately get the result 

{deti[(Fs)cot(Fs)]}'/' = 

= gcoth (^s^J'+ \f¥'^^Q'^ cot [s^-T^^T'^^Q'^ 



(162) 
(163) 



(164) 



We now evaluate {(Fs) cot(Fs)}oo- We define 
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and a sign 



Then, the matrices 



s:=QI\Q\ = \Q\IQ 



A — ^ in 



2(// + //) 
satisfy the relations 

Fii\A(^i)Xu = fi^{i)ijiv = ^{i),j.xF\v and ^ A(^i)ixv = 

i 

Using these two relations, we have for a function 



(165) 
(166) 

(167) 

(168) 



This is the key point of introducing the matrices Now, from this relation, we have 

{{Fs) cot(eFs)}oo = HMfis) cot(/is) . (169) 



Since 



we have 



and thus 



Foo = = Foo and FqxFxq = E^ 



(170) 
(171) 



{(F5)cot(Fs)}oo = 7^372 {(/- - E^)f +cot{Us) - {fl - i2)/_cot(/_s)} . (172) 



We have 



and 



We also have 



(173) 
(174) 



/+cot(/+s) = V ^ + + ^2 coth(s V ^ + + ^2) 
/_ cot(/_s) =\J-J^+^'^ + g^ coi{s\l -J^+^J'^ + g^) 



(175) 



So we get 



{eFcot(eFs)}oo 



I - ^(52 + ^2) + + g2|g^jr + + gicoth (^es-\/j^+ V^2 + g2^ 



(176) 
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C Treatment of Large Chemical Potential 



In Section r2.5.H we have pointed out that a special treatment is necessary when the chemical 
potential of the system is too large. In this appendix, we give the details of how to handle 
this case. The discussion here is somewhat similar to Reference [l2j. For simplicity, we first 
consider the system without the constant background field, hence the goal is to evaluate 



4vr3/2 



dss ^ exp [ 



s{{27tTI -ifi)'^ + m'^}] , 



(177) 



for a general value of The inclusion of -B-field (and i^-field) can be done in relative ease. 

We first discuss the analytic continuation of the propertime in general. Recall that we 
have dealt with the "Euclidean worldline Lagrangian" appearing in Equations (I33|34p . This 
stems from our choice of the propertime variable in Equation (j29p . There, we chose it so 
that the exponential appears as real, whereas Schwinger chooses in Ref. [15] as 



r-oo —iXt 
\nX = - I —^dt . 



(178) 



This choice leads to the evolution equation of the real propertime. Therefore, we can think 
of our choice as the analytically continued version, i.e., the imaginary propertime s. See 
Figure [ 



\ 



Figure 39: Analytic continuation of real to imaginary propertime. Since our notation is that s is the 
imaginary propertime, the vertical axis corresponds to the real propertime. 

We note that the situation is a little more complicated when we must introduce the 
cutoff at s = 1. In this case, it seems as if we have discrepancy with the real propertime 
case due to the quarter circle contribution. See Figure HOI However, the cutoff is designed 
to make the imaginary propertime integral well-defined, and it is nonsensical in the real 
one. Therefore, we define the world of real propertime as the analytic continuation of 
the imaginary propertime with the cutoff. This implies that the contour Ci in Figure HOl 
corresponds to Cr + Cps. 

Now, we notice that we have choices to analytically continue from the real to imaginary 
propertime: t ^ s or t — > — s, in other words, it is also possible to analytically continue to 
the left of Figure [39l The former choice has been made in Equation (I29p and we found that 
this choice is not suitable when the chemical potential is too large. When this is the case. 
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Cr 




C1 



Figure 40: With the introduction of the cutofT, the real and imaginary propertimes have difference due 
to the quarter circle. We therefore, define the real counterpart with the quarter-circle contribution: Ci <4> 

Cr + Cp3 . 



I 




Figure 41: The analytic continuation from positive to negative real axis. The contour Ci now corresponds 
to C2 + C3, where C3 := Cp3 + CmS- 

we must pick the other choice. This amounts to the analytic continuation of the variable s 
in Equation (|177|) from the positive to the negative direction as shown in Figure BTl 

We must be careful about which way we should close the contour. The imaginary part 
of the exponential in Equation p77|) is 2Tii^Tls and this tells us that we should close the 
contour above when / is positive and below when I is negative, so that the large circles do 
not contribute to the integral. See Figure W2\ Finally, we must be careful with the branch 
cut that arise from the factor of s~^/^ in Equation (|177p . The contour should not cross the 
cut and this is possible only when the cut runs along the negative horizontal axis. This 
causes the ambiguity in choosing the branch of on the negative real axis. We take 

the convention that when we close the contour below, we choose the second branch value. 
Later, we will show that this is the consistent choice. 

When we have the B field, we have poles along the vertical axis in the s-plane and we 
must pick them up, as they are within the contour. We will deal with this later. 

Let us actually evaluate Equation (11770 . As mentioned in Section 12. 5. H the integral 
converges only when 

{2'KTlf - /u^ + > , (179) 
is satisfied. Without loss of generality, let us assume that fi > and m > 0. When fi < m, 
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Figure 42: Analytic continuation of the variable s from positive to negative real axis. The upper (resp. 
lower) contour should be chosen for positive (resp. negative) I. The branch cut comes from the factor of 
s~^^^ in Equation (I177p . The branch cut must run along the negative real axis to avoid the cross-over with 
the contours. Our convention is that when we take the lower contour, the value of s"^''^ is taken to be on 
the second branch for the negative real axis. 



this condition is satisfied for all I, but when /i > m, it is satisfied only for |Z| > / with 



/ := [{2TTT)-^e{fi - m)^y fi"^ - m^]^ , (180) 

where 0{x) is the Heaviside theta and the symbol [x]g is the half-integer version of the Gauss 
symbol, i.e., it is the largest half integer less than or equal to x. Notice, in particular, that 
when fi < m, we have I = —1/2 and / runs over all half integers. According to the previous 
discussion. Equation p77|) must be split at / and for positive and negative /. Thus, we write 



Eqn. (fTTTD = I dss-^/"^ ^ exp [ - s{(27rr/ - i^if + n?]] 

■'^^ \i\>i 

+ [ dss-^'^ y exp [ - s{{2t:TI - ifif + mH] 

T r 

+ -^/ ^ exp [-s{(27rr/-i^)2+m2}] , (181) 

47r J (C9+C:\]down , r 



where "up" and "down" indicated on the integration contour implies the way the contour 
is closed as in Figure H2l Consider the integral in the last line. Since the both integral and 
sum are well-defined now, we may exchange them. Also we send the summation variable I 
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to — /. We then have 



4vr3/2 



(=1/2 



(C2+C3)down 



dsfi-^/^exp [ - s{{2TTTl + iijf + m^}] 



T 



, fCsMown 



4vr3/2 
T 



«=l/2 



A 3/2 Yj dss^^^/^exp[s{{27rTl + inf + m^}] 



. r ^ 

Z=l/2 ^ 



iei*'^ exp [ - e-*<^{(27rT/ + ifif + m^}] 



(182) 



where we have taken the second branch value for the negative real axis as our convention 
and the sign of the integral variable (p was reversed so that now the integration range is (0, it) 
rather than (0, — vr). A little inspection shows that this is exactly the complex conjugate of 
the second line in Equation (jlSlh . 

We can also split the first line in Equation (jl8ip into positive and negative sums and it 
is easy to see that they are complex conjugate of each other. Therefore, we can write 



Eqn. (fTTTD 



27r3/2 
I 



1>T •'^ 

/oo 
dss-^'^ exp [s{(27rr/ - i^if + m^}] 

+ [ #e-i*'^exp[-e*'^{(2^rZ-i/x)2 + m2}] 



1=1/2 ■ 



(183) 



After so much ado, we point out the following. In Equation ()177p . assuming that |/| > /, 
we can, in principle, carry out s-integral before the sum. The result is a function of / and 
valid only when |/| > I. However, we can analytically continue the integral beyond the valid- 
ity range of / to cover all half integers, without worrying about the case by case examination 
that we have done. In fact, one can find a closed form of the s integral in Equation (I177p 
and verify this claim by comparing with the result (|183p . This in turn implies the validity 
of the choice of the branch cut for the lower contour discussed before. In practice, when 
the B field is present, we are not able to find the closed form of s-integral and must resort 
to the case by case integral. In what follows, we discuss this case more. 

As mentioned earlier, the case with the pure -B-field has simple poles along the imaginary 
axis in the s-plane, due to the factor of hyperbolic cosine in Equation ([571) . Therefore, 
when we must analytically continue s to the negative horizontal axis, we need to pick up 
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the poles. The condition for I, which determines the point for the analytic continuation, 
is not affected by the presence of the hyperbolic cosine, because it approaches 1 at large 
s. Hence the expressions for this case is similar to the previous case except for the extra 
contributions from the poles. 

Zs = 

T ^ 
1 Re 



2r 



+ i 



ds{sB) coth(sS)s-^/2 ^ gxp [ - s{(27rT/ - i\if + m^}] 
l>l 

/oo ' 
ds{sB) coth(s5)s-^/2 ^ g^p [s{(27rT/ - i^if + m^}] 

/=l/2 

r d(\>{e''^B) coth(e^'^S)e-i^'^ V exp [ - e^'^{(27rrZ - i^if + m^}] 

-'0 , wo 



/=l/2 



i oo 



l=l/2k=ki ^ ^ 

Since the location of the poles are at 

s = inir/B for all n G Z , 



(184) 



(185) 



for sufficiently large S, the pole at n = 1 (and more) comes inside the unit circle that we 
are excluding. Thus naively, we might exclude those poles that are inside the unit circle, 
i.e., we set 



B 



TT 



+ 1 



(186) 



J G 



where [yjc here is the usual integer Gauss symbol. However, numerical evaluation reveals 
that this cross over results in a huge jump in the effective potential and it does not appear 
physical. This is the indication that the B field may not be much larger than the cutoff so 
that the poles of n > should not cross the circle. This means that we must have B < tt 
and kb = 1. 
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